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Abstract. It is widely believed that the celebrated 2D Ising model at criticality 
has a universal and conformally invariant scaling limit, which is used in deriving 
many of its properties. However, no mathematical proof has ever been given, and 
even physics arguments support (a priori weaker) Mobius invariance. We introduce 
discrete holomorphic fermions for the 2D Ising model at criticality on a large family 
of planar graphs. We show that on bounded domains with appropriate boundary 
conditions, those have universal and conformally invariant scaling limits, thus proving 
the universality and conformal invariance conjectures. 
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1. Introduction 
1.1. Universality and conformal invariance in the Ising model. 

1.1.1. Historical background. The celebrated Lenz-Ising model is one of the simplest 
systems exhibiting an order-disorder transition. It was introduced by Lenz in |Len20j . 
and his student Ising proved [Isi25j in his PhD thesis the absence of phase transition 
in dimension one, wrongly conjecturing the same picture in higher dimensions. This 
belief was widely shared, and motivated Heisenberg to introduce his model |Hei28] . 
However, some years later Peierls |Pei36] used estimates on the length of interfaces 
between spin clusters to disprove the conjecture, showing a phase transition in the two 
dimensional case. After Kramers and Wannier |KW41] derived the value of the critical 
temperature and Onsager |Ons44j analyzed behavior of the partition function for the 
Ising model on the two-dimensional square lattice, it became an archetypical example 
of the phase transition in lattice models and in statistical mechanics in general, see 
[NisOSl INisOQj for the history of its rise to prominence. 

Over the last six decades, thousands of papers were written about the Ising model, 
with most of the literature, including this paper, restricted to the two dimensional case 
(similar behavior is expected in three dimensions, but for now the complete description 
remains out of reach). The partition function and other parameters were computed 
exactly in several different ways, usually on the square lattice or other regular graphs. It 
is thus customary to say that the 2D Ising model is exactly solvable, though one should 
remark that most of the derivations are non-rigorous, and moreover many quantities 
cannot be derived by traditional methods. 

Arrival of the renormalization group formalism (see |Fis98j for a historical ex- 
position) led to an even better physical understanding, albeit still non-rigorous. It 
suggests that block-spin renormalization transformation (coarse-graining, i.e., replac- 
ing a block of neighboring sites by one) corresponds to appropriately changing the scale 
and the temperature. The Kramers- Wannier critical point arises then as a fixed point 
of the renormalization transformations, with the usual picture of stable and unstable 
directions. 

In particular, under simple rescaling the Ising model at the critical temperature 
should converge to a scaling limit - a "continuous" version of the originally discrete 
Ising model, which corresponds to a quantum field theory. This leads to the idea of 
universality, the Ising models on different regular lattices or even more general planar 
graphs belong to the same renormalization space, with a unique critical point, and so 
at criticality the scaling limit and the scaling dimensions of the Ising model should be 
independent of the lattice (while the critical temperature depends on it). Being unique, 
the scaling limit at the critical point is translation and scale invariant, which allows to 
deduce some information about correlations [PP66[ [Kad 66]. By additionally postulat- 
ing invariance under inversions, one obtains Mobius invariance, i.e. invariance under 
global conformal transformations of the plane, which allows |Pol70] to deduce more. 
In seminal papers [BPZ84a| IBPZ84b] Belavin, Polyakov and Zamolodchikov suggested 
much stronger full conformal invariance (under all conformal transformations of subre- 
gions), thus generating an explosion of activity in conformal field theory, which allowed 
to explain non-rigorously many phenomena, see |ISZ88] for a collection of the founding 
papers of the subject. Note that in the physics literature there is sometimes confusion 
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between the two notions, with Mobius invariance often called conformal invariance, 
though the latter is a much stronger property. 

Over the last 25 years our physical understanding of the 2D critical lattice mod- 
els has greatly improved, and the universality and conformal invariance are widely 
accepted by the physics community. However, supporting arguments are largely non- 
rigorous and some even lack physical motivation. This is especially awkward in the 
case of the Ising model, which indeed admits many exact calculations. 

1.1.2. Our results. The goal of this paper is to construct lattice holomorphic fermions 
and to show that they have a universal conformally invariant scaling limit. We give 
unambiguous (and mathematically rigorous) arguments for the existence of the scaling 
limit, its universality and conformal invariance for some observables for the 2D Ising 
model at criticality, and provide the framework to establish the same for all observ- 
ables. By conformal invariance we mean not the Mobius invariance, but rather the full 
conformal invariance, or invariance under conformal transformations of subregions of 
C. This is a much stronger property, since conformal transformations form an infinite 
dimensional pseudogroup, unlike the Mobius ones. Working in subregions necessarily 
leads us to consider the Ising model in domains with appropriate boundary conditions. 

At present we cannot make rigorous the renormalization approach, but we hope 
that the knowledge gained will help to do this in the future. Rather, we use the 
integrable structure to construct discrete holomorphic fermions in the Ising model. 
For simplicity we work with discrete holomorphic functions, defined e.g. on the graph 
edges, which when multiplied by the \fdz field become fermions or spinors. Those 
functions turn out to be discrete holomorphic solutions of a discrete version of the 
Riemann-Hilbert boundary value problem, and we develop appropriate tools to show 
that they converge to their continuous counterparts, much as Courant, Friedrichs and 
Lewy have done in |CFL28] for the Dirichlet problem. The continuous versions of our 
boundary value problems are \/d^-covariant, and conformal invariance and universality 
then follow, since different discrete conformal structures converge to the same universal 
limit. 

Starting from these observables, one can construct new ones, describe interfaces 
by the Schramm's SLE curves, and prove and improve many predictions originating 
in physics. Moreover, our techniques work off criticality, and lead to massive field 
theories and SLEs. Several possible developments will be the subject of our future 
work [CSm ISmilli iKSTTl HCT] . 

We will work with the family of isoradial graphs or equivalently rhombic lattices. 
The latter were introduced by Duffin |Duf68] in late sixties as (perhaps) the largest 
family of graphs for which the Cauchy-Riemann operator admits a nice discretization. 
They reappeared recently in the work of Mercat [MerOlj and Kenyon |Ken02j , as iso- 
radial graphs - possibly the largest family of graphs were the Ising and dimer models 
enjoy the same integrability properties as on the square lattice: in particular, the crit- 
ical point is well defined, with weights depending only on the local structure. More 
recently, Boutilier and de Tiliere |BdT08t IBdTOQ] used the Fisher representation of 
the Ising model by dimers and Kenyon's techniques to calculate, among other things, 
free energy for the Ising model on isoradial graphs. While their work is closely related 
to ours (we can too use the Fisher representation instead of the vertex operators to 
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construct holomorphic fermions), they work in the full plane and so do not address 
conformal invariance. Note that earlier eight vertex and Ising models were considered 
by Baxter |Bax78] on Z-invariant graphs, arising from planar line arrangements. Those 
graphs are topologically the same as the isoradial graphs, and the choice of weights 
coincides with ours, so quantities like partition function would coincide. Kenyon and 
Schlenker |KS05j have shown that such graphs admit isoradial embeddings, but those 
change the conformal structure, and one does not expect conformal invariance for the 
Ising model on general Z-invariant graphs. 

So there are two reasons for our choice of this particular family: firstly it seems to 
be the largest family where the Ising model we are about to study is nicely defined, and 
secondly (and perhaps not coincidentally) it seems to be the largest family of graphs 
where our main tools, the discrete complex analysis, works well. It is thus natural to 
consider this family of graphs in the context of conformal invariance and universality 
of the 2D Ising model scaling limits. 

The fermion we construct for the random cluster representation of the Ising model 
on domains with two marked boundary points is roughly speaking given by the proba- 
bility that the interface joining those points passes through a given edge, corrected by 
a complex weight. The fermion was proposed in |Smi06[ ISmiOTj for the square lattice 
(see also independent |RC06] for its physical connections, albeit without discussion 
of the boundary problem and covariance). The fermion for the spin representation is 
somewhat more difficult to construct, it corresponds to the partition function of the 
Ising model with a ^/z monodromy at a given edge, again corrected by a complex 
weight. We describe it in terms of interfaces, but alternatively one can use a product 
of order and disorder operators at neighboring site and dual site, or work with the 
inverse Kasteleyn's matrix for the Fisher's dimer representation. It was introduced in 
[SmiOGj . although similar objects appeared earlier in Kadanoff and Ceva [KC71j (with- 
out complex weight and boundary problem discussions) and in Mercat [MerOlj (again 
without discussion of boundary problem and covariance). 

Complex analysis on isoradial graphs is more complicated then on the square grid, 
and less is known a priori about the Ising model there. As a result parts of our paper 
are quite technical, so we would recommend reading the much easier square lattice 
proofs |Smi07l ICSllj . as well as the general exposition |Smi06[ ISmilOj first. 

1.1.3. Other lattice models. Over the last decade, conformal invariance of the scaling 
limit was established for a number of critical lattice models. An up-to-date introduction 
can be found in |Smi06 j. so we will only touch the question of universality here. 

Spectacular results of Kenyon on conformal invariance of the dimer model, see e.g. 
|KenOOa[ IKenOlj . were originally obtained on the square lattice. Some were extended 
to the isoradial case by de Tihere |dT07] . but the questions of boundary conditions 
and hence conformal invariance were not addressed yet. 

Kenyon's dimer results had corollaries [KenOObj for the Uniform Spanning Tree 
(and the Loop Erased Random Walk). Those used the Temperley bijection between 
dimer and tree configurations on two coupled graphs, so they would extend to the 
situations where boundary conditions can be addressed and Temperley bijection exists. 

Lawler, Schramm and Werner used in |LSW04j simpler observables to establish 
conformal invariance of the scaling limit of the UST interfaces and the LERW curves. 
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and to identify them with Schramm's SLE curves. In both cases one can obtain observ- 
ables using the Random Walk, and for the UST one can use the Kirchhoff circuit laws 
to obtain discrete holomorphic quantities. The original paper deals with the square 
lattice only, but it easily generalizes whenever boundary conditions can be addressed. 

In all those cases we have to deal with convergence of solutions of the Dirichlet 
or Neumann boundary value problems to their continuous counterparts. While this is 
a standard topic on regular lattices, there are technical difficulties on general graphs; 
moreover functions are unbounded (e.g. observable for the LERW is given by the 
Poisson kernel), so controlling their norm is far from trivial. Tools developed by us in 
jCS08] for use in the current paper however resolve most of such difficulties. 

Situation is somewhat easier with the observables for the Harmonic Explorer and 
Discrete Gaussian Free Field, as discussed by Schramm and Sheffield |SS05| ISS09] - 
both are harmonic and solving Dirichlet problem with bounded boundary values, so 
the generalization from the original triangular lattice is straight-forward. Note though 
that the key difficulty in the DGFF case is to establish the martingale property of the 
observable. 

Unlike the observables above, the one used for percolation in |Smi01a[ ISmiOlb] is 
very specific to the triangular lattice, so the question of universality is far from being 
resolved. 

All the observables introduced so far (except for the fermions from this paper) are 
essentially bosonic, either invariant under conformal transformations (p or changing like 
"pre-pre-Schwarzian" forms, i.e. by an addition of const ■ (p' . They all satisfy Dirichlet 
or Neumann boundary conditions, when establishing convergence is a classical subject, 
dating back to Courant, Friedrichs and Lewy |CFL28j . albeit in the non-bounded case 
one meets serious difficulties. 

In the Ising case we work with fermions, hence the Riemann-Hilbert boundary 
value problem (or rather its homogeneous version due to Riemann). Such problems turn 
out to be much more complicated already on regular lattices: near rough boundaries 
(which arise naturally since interfaces are fractal) our observables blow up fast. When 
working on general graphs, the main problem remains, but the tools become quite 
limited. 

We believe that further progress in other models requires the study of holomorphic 
parafermions [SmiOGj . so we expect even more need to address the Riemann boundary 
value problems in the future. 

1.2. Setup and main results. Throughout the paper, we work with isoradial graphs 
or, equivalently, rhombic lattices. A planar graph F embedded in C is called 6-isoradial 
if each face is inscribed into a circle of a common radius 6. If all circle centers are 
inside the corresponding faces, then one can naturally embed the dual graph F* in 
C isoradially with the same 6, taking the circle centers as vertices of F*. The name 
rhombic lattice is due to the fact that all quadrilateral faces of the corresponding 
bipartite graph A (having F U F* as vertices and radii of the circles as edges) are 
rhombi with sides of length 6. We denote the set of rhombi centers by (} (example of 
an isoradial graph is drawn in Fig. [IJA.). We also require the following mild assumption: 

the rhombi angles are uniformly bounded away from and vr 
(in other words, all these angles belong to [r],TT — r]] for some fixed ?7 > 0). Below we 
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(A) (C) 

Figure 1. (A) Example of an isoradial graph F (black vertices, solid 
lines), its dual isoradial graph F* (gray vertices, dashed lines), the corre- 
sponding rhombic lattice or quad-graph (vertices A = F U F*, thin lines) 
and the set (} = A* (rhombi centers, white diamond-shaped vertices). 
(B) Local notation near m G F, with neighbors of u enumerated coun- 
terclockwise by 1, 2 . . . , s, s -|- 1, . . . , n. The weight /xf(w) is equal to the 
shaded polygon area. (C) Definition of s-holomorphic functions: F{zq) 
and F{zi) have the same projections on the direction [i{w—u)]^2. Thus, 
we have one real identity for each pair of neighboring zq, z^. 

often use the notation const for absolute positive constants that don't depend on the 
mesh b or the graph structure but, in principle, may depend on r^. We also use the 
notation f ^ g which means that a double-sided estimate consti -f^g^ const2 -g 
holds true for some consti^2 > which are independent of S. 

It is known that one can define the (critical) Ising model on F* so that 

(a) the interaction constants J^-^^ are local (namely, depend on the lengths of edges 
connecting Wij G F*, ~ Wj, only) and 

(b) the model is invariant under the star-triangle transform. 

Such invariance is widely recognized as the crucial sign of the integrability. Note 
that the star-triangle transform preserves the isoradial graph/rhombic lattice structure. 
Moreover, isoradial graphs form the largest family of planar graphs (embedded into C) 
satisfying these properties (see [CS06] and references therein). At the same time. 
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discrete holomorphic functions on isoradial graphs provide the simplest example of 
a discrete integrable system in the so-called "consistency approach" to the (discrete) 
integrable systems theory (see |BS08] ). 

Recently, it was understood (see |Mer01t ISmiOGl [RC06] ) that some objects coming 
from the theoretical physics approach to the Ising model (namely, products of order and 
disorder operators with appropriate complex weights) can be considered and dealt with 
as discrete holomorphic functions (see Sect. l3.2l for further discussion). These functions 
(which we call basic observables or holomorphic fermions) provide a powerful tool for 
rigorous proofs of several results concerning the conformal invariance of the critical 
2D Ising model. Implementing the program proposed and started in |Smi06| ISmi07] 
for the square grid, in this paper we mainly focus our attention on the topologically 
simplest case, when the model is defined in the simply-connected (discrete) domain 
having two marked boundary points a^, (but see Section [6] for more involved setup). 
We have to mark some boundary points so that the conformal modulus is non-trivial, 
allowing us to construct conformal invariants. 

We will work with two representations of the Ising model: the usual spin, as well 
as the random cluster (Fortuin-Kasteleyn). The observables are similar, but do not 
directly follow from each other, and require slightly different approaches. In either case 
there is an interface (between spin clusters or random clusters) - a discrete curve 7*^ 
running from to b^ inside (see Section [2m 12.21 for precise definitions). In both 
cases the basic observables are martingales with respect to (filtration induced by) the 
interface grown progressively from a^, which opens the way to identify its scaling limit 
as a Schramm's SLE curve, cf. |Smi06] . 

The interface in the random cluster representation can in principle pass through 
some point twice, but with our setup we move apart those passages, so that the curve 
becomes simple and when arriving at the intersection it is always clear how to proceed. 
This setup is unique where the martingale property holds, so there is only one con- 
formally invariant way to address this problem. Note that the resulting scaling limit, 
the SLE(16/3) curve, will have double points. A similar ambiguity arises in the spin 

model (when, e.g. a vertex is surrounded by four spins " — I h"), but regardless of 

the way to address it (e.g. deterministic, like always turning right, or probabilistic, like 
tossing a coin every time) the martingale property always holds, and so the SLE(3) is 
the scaling limit. The latter is almost surely simple, so we conclude that the double 
points in the discrete case produce only very small loops, disappearing in the scaling 
limit. 

The first two results of our paper say that, in both representations, the holomor- 
phic fermions are uniformly close to their continuous conformally invariant counter- 
parts, independently of the structure of (or <)'^) and the shape of fl^ (in particular, 
we don't use any smoothness assumptions concerning the boundary). Namely, we prove 



the following two theorems, formulated in detail as Theorems 14.31 and 15.61 



Theorem A (FK-Ising fermion). Let discrete domains {n^;a^,b^) with two marked 
boundary points a^,b^ approximate some continuous domain {Q;a,b) as 6 0. Then, 
uniformly on compact subsets of Q and independently of the structure of (}^ , 




UNIVERSALITY IN THE 2D ISING MODEL 



9 



where F^[z) = F^{z;Q^ ,b^) is the discrete holomorphic fermion and $ denotes the 
conformal mapping from Q onto the strip M x (0, 1) such that a, b are mapped to =Foo. 



Theorem B (spin-Ising fermion). Let discrete domains {fl^;a^,b^) approximate some 
continuous domain {Q;a,b) as 5 — )■ 0. Then, uniformly on compact subsets of Q and 
independently of the structure ofT^, 



where = F\z; Q\ a\ b^) is the discrete holomorphic fermion and \1/ : f2 — t- C+ is the 
conformal mapping such that a and b are mapped to oo and 0, appropriately normalized 
at b. 

Because of the aforementioned martingale property, these resuhs are sufficient to 
prove the convergence of interfaces to conformally invariant Schramm 's SLE curves (in 
our case, SLE(3) for the spin representation and SLE(16/3) for the FK representation) 
in the weak topology given by the convergence of driving forces in the Loewner equation, 
cf. |LSW04[ ISmiOGj . A priori this is very far from establishing convergence of curve 
themselves, but we use techniques of |KS09j to prove this stronger convergence in 



The third result shows how our techniques can be used to find the (conformally 
invariant) limit of some macroscopic quantities, "staying on the discrete level", i.e. 
without consideration of the limiting curves. Namely, we prove a crossing probability 
formula for the critical FK-Ising model on isoradial graphs, analogous to Cardy's for- 
mula [SmiOlal ISmiOlbj for critical percolation and formulated in detail as Theorem 16. II 

Theorem C (FK-Ising crossing probability). Let discrete domains {Q^]a^,b^,c^,d^) 
with alternating (wired/ free/ wired/ free) boundary conditions on four sides approximate 
some continuous topological quadrilateral {Q; a, b,c,d) as 6 ^ 0. Then the probability of 
an FK cluster crossing between two wired sides has a scaling limit, which depends only 
on the conformal modulus of the limiting quadrilateral, and is given for the half-plane 
by 



The version of this formula for multiple SLEs was derived by Bauer, Bernard and 
Kytola in |BBK05j . see page 1160, their notation for the modulus related to ours by 
X = 1 — u. Besides being of an independent interest, this result together with [KS09j 
is needed to improve the topology of convergence of FK-Ising interfaces. Curiously, 
the (macroscopic) answer for a unit disc (D; — e*'^, e"*"^, e**^, —e~^'^) formally coincides 
with the relative weights corresponding to two possible crossings inside (microscopic) 
rhombi (see Fig. |2JA.) in the critical model (see Remark I6.2p . 

1.3. Organization of the paper. We begin with the definition of Fortuin-Kasteleyn 
(random cluster) and spin representations of the critical Ising model on isoradial graphs 
in Section |2j From the outset we work with critical interactions, but in principle one 
can introduce a temperature parameter, which would lead to massive holomorphic 
fermions. We also introduce the basic discrete holomorphic observables (holomorphic 



and 




(1.1) 
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fermions) satisfying the martingale property with respect to the growing interface and, 
essentially, show that they satisfy discrete version of the Cauchy-Riemann equation 
(Proposition 12.21 and Proposition 12. 5p using some simple combinatorial bijections be- 
tween the sets of configurations. Actually, we show that our observables satisfy the 
stronger "two-points" equation which we call spin or strong holomorphicity, or simply 
s-holomorphicity. 

We discuss the properties of s-holomorphic functions in Section [31 The main 
results are: 

(a) The (rather miraculous) possibility to define naturally the discrete version of 
h{z) = Im J{f{z)ydz, see Proposition 13.61 Note that the square (/(z))^ of a 
discrete holomorphic function f{z) is not discrete holomorphic anymore, but 
unexpectedly it satisfies "half" of the Cauchy-Riemann equations, making its 
imaginary part a closed form with a well-defined integral; 

(b) The sub- and super-harmonicity of h on the original isoradial graph F and its 
dual r*, respectively, and the a priori comparability of the components 
and /i|p, which allows one to deal with h as with a harmonic function: e.g. 
nonnegative /I's satisfy a version of the Harnack Lemma (see Section [33]); 

(c) The uniform (w.r.t. 6 and the structure of the isoradial graph/rhombic lattice 
r,r*/<)) boundedness and, moreover, uniform Lipschitzness of s-holomorphic 
functions inside their domains of definition , with the constants depending 
on M = maXj^gQ* \h{v)\ and the distance d = dist(2;; dQ^) only (Theorem 13. 12^ 
these results should be considered as discrete analogous of the standard esti- 
mates from the classical complex analysis); 

(d) The combinatorial trick (see Section 13. 6p that allows us to transform the dis- 
crete version of the Riemann-type boundary condition f{() \\ {t{())~2 into the 
Dirichlet condition for h^g^ on both F and F*, thus completely avoiding the ref- 
erence to Onsager's magnetization estimate used in |Smi06[ ISmiOTj to control 
the difference — on the boundary. 

We prove the (uniform) convergence of the basic observable in the FK-Ising model 
to its continuous counterpart in Section [H The main result here is Theorem 14. 3[ which 
is the technically simplest of our main theorems, so the reader should consider the proof 
as a basic example of the application of our techniques. Besides the results from [CS08j 
and the previous Sections, the important idea (exactly as in |Smi07j ) is to use some 
compactness arguments (in the set of all simply-connected domains equipped with the 
Caratheodory topology) in order to derive the uniform (w.r.t. to the shape of fl^ and 
the structure of <C>'^) convergence from the "pointwise" one. 

In Section [5] we prove analogous convergence result for the holomorphic fermion 
defined for the spin representation of the critical Ising model (Theorem 15. 6p . There are 
two differences from the preceding Section: the unboundedness of the (discrete) integral 
h = Im J {f{z)Ydz (this prevents us from the immediate use of compactness arguments) 
and the need to consider the normalization of our observable at the target point 
(this is crucial for the martingale property). In order to handle the normalization at 
6'', we assume that our domains VL^ contain a (macroscopic) rectangle near and their 
boundaries dVL^ approximate the corresponding straight segment as 5 — )■ 0. Making 
this technical assumption, we don't lose much generality, since the growing interface. 
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though fractal in the hmit, doesn't change the shape of the domain near b^. Then, we 
use a version of the boundary Harnack principle (Proposition I5.3p in order to control 
the values of h in the bulk through the fixed value f{h^)- Another important technical 
ingredient is the universal (w.r.t. to the structure of <)'') multiplicative normalization 
of our observable. Loosely speaking, we define it using the value at of the discrete 
holomorphic fermion in the discrete half-plane (see Theorem 15.41 for further details). 

Section E] is devoted to the crossing probability formula for the FK-Ising model on 
discrete quadrilaterals [VL^] a*^, 6'', c'', d^) (Theorem 16. 11 see also Remark l6.2l) . The main 
idea here is to construct some discrete holomorphic in VL^ function whose boundary 
values reflect the conformal modulus of the quadrilateral. Namely, in our construction, 
discrete functions = Im J [f^[z)Ydz approximate the imaginary part of the confor- 
mal mapping from VL^ onto the slit strip [R x (0, 1)] \ (— cxo + iK\ ix] such that is 
mapped to the "lower" —oo, to -|-oo; to the "upper" —oo and to the tip i>i. The 
respective crossing probabilities are in the 1-to-l correspondence with values which 
approximate x as 5 — t- 0. Since x is uniquely determined by the limit of conformal 
moduli of {Q^; a^, b\ c\ d^), we obtain ([Ll]) (see further details in Section [6]). Finally, 
Appendix contains several auxiliary lemmas: estimates of the discrete harmonic mea- 
sure, discrete version of the Cauchy formula, and technical estimates of the Green 
function in the disc. We refer the reader interested in a more detailed presentation of 
the discrete complex analysis on isoradial graphs to our paper |CS08j . 
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2. Critical spin- and FK-Ising models on isoradial graphs. 
Basic observables (holomorphic fermions) 

2.L Critical FK-Ising model. 

2.LL Loop representation of the model, holomorphic fermion, martingale property. 
We will work with a graph domain which can be thought of as a discretization of a 
simply-connected planar domain with two marked boundary points. Let C <C> be a 
simply-connected discrete domain composed of inner rhombi z G Int and boundary 
half-rhombi ( G dQ^^, with two marked boundary points a^, b^ and Dobrushin boundary 
conditions (see Fig. [2]B): dQ^^ consists of the "white" arc a^^t; ^^e "black" arc b^a^, 
and two edges [a^a^], [&b^t] of A. Without loss of generality, we assume that 

6^ — 6^ = i6, i.e., the edge b^ = ^■^ oriented vertically. 

For each inner rhombus z G Int fi^ we choose one of two possibilities to connect its 
sides (see Fig. [2J\, there is only one choice for boundary half-rhombi), thus obtaining 
the set of configurations (whose cardinality is 2*^^'^*'^0''). The partition function of 



12 



DMITRY CHELKAK AND STANISLAV SMIRNOV 




(A) (B) 

Figure 2. (A) Loop representation of the critical FK-Ising model 
on isoradial graphs: the relative weights corresponding to two possible 
choices of connections inside the inner rhombus z (the partition function 
is given by fl2.2p ). (B) Discrete domain with a sample configuration. 
Besides loops, there is an interface 7'' connecting to . Calculating the 
winding(7''; 6'^ w we draw 7^ so that it intersects the edge ^ = [^b'Cw] 
orthogonally. As 7'' grows, it separates some part of (shaded) from 
h^. We denote by ^% \ [a^7f-7j] 

the connected component containing h 

(unshaded) . 



the critical FK-Ising model is given by 

Z=Y^ ^2*^^°°^^^ n.eintf.*, sin i ^(z), (2.2) 

config. 

where 6'config. (2;) is equal to either 6 ot 6* = ^ — depending on the choice of connections 
inside rhombus z (see Fig. [2]A.). 

We described the loop representation, since at criticality it is easier to work with, 
than the usual random cluster one. The loops trace the perimeters of random clusters, 
and the curve joining the two marked boundary points is the interface between a 
cluster and a dual cluster wired on two opposite boundary arcs (the so-called Dobrushin 
boundary conditions). 

Let ^ = [^b^w] be some inner edge of (where ^ T, ^ T*)- Due to the 
boundary conditions chosen, each configuration consists of (a number of) loops and 
one interface 7*^ running from to b^. The holomorphic fermion is defined as 



^'(0 = i^(Ua^6*)(0:=(25)-^-E ;^(ee7^)-e-^-°^-s(^'''"-5^ , (2.3) 
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where 

the indicator function of the event that the interface intersects ^ and 

winding(7''; 6'^ -w ^) = winding(7'^; a'' -w ^) — winding(7'^; -w b^) (2.4) 

denotes the total turn of 7'' measured (in radians) from to ^. Note that, for all 
configurations and edges ^ one has (see Fig. [2]B), 

Remark 2.1 (Martingale property). For each^, F^q^i,^^^/,^^ ^^sy^s o martin- 

gale with respect to the growing interface [a^ = 7o,7i, •••57^5 •••) (tM the stopping time 
when 7*^ hits ^ or ^ becomes separated from b^ by the interface, see Fig. \^). 

Proof. Since the winding(7''; b^ ^) doesn't depend on the beginning of the interface, 
the claim immediately follows from the total probability formula. □ 



2.1.2. Discrete boundary value problem for F^ . We start with the extension of F^ to 



the centers of rhombi z G fil. Actually, the (rather fortunate) opportunity to use 



the definition given below reflects the discrete holomorphicity of F^ (see discussion in 
Section [3I2I) . 

Proposition 2.2. Let z G Int be the center of some inner rhombus UiWiU2W2- 
Then, there exists a complex number F^{z) such that 

F\[u,Wk]) = PToi[F\z); [t{wk-u,)]-"^^ , 3,k = 1,2. (2.5) 

The proposition essentially states that F^ is spin holomorphic as specified in Defini- 
tion [XT] below. By Proj[X; u] we denote the orthogonal projection of the vector X on 
the vector u, which is parallel to u and equal to 



Proj X; M = Re X— — = - + 
\ 1^1/ 1^1 2 



2\u\ 



(here we consider complex numbers as vectors). Because of the latter rewriting, the 
choice of the sign in the square root in (12.51) does not matter. 

Proof. As on the square grid (see |Smi07j ). the proof is based on the bijection between 
configurations which is produced by their local rearrangement at z. It is sufficient 
to check that the contributions to F^ {[ujW klYs of each pair of configurations drawn 
in Fig. [3] are the specified projections of the same complex number. The relative 
contributions of configurations (up to the same real factor coming from the structure 
of the configuration away from z) to the values of F on four edges around z for the 
pairs "L" and "R" are given by 





e'2. F^{[u2Wi]) 


e'2 . F\[U2W2]) 


e*2 ■ F^{[uiW2]) 


e*2 ■ F^{[uiWi]) 


L 


sin 1 + sinf 


e*^ sin f 


—i sin f 


e-^^*[72sinf + sinf ] 


R 


sinf + 72 sinf 


e^^[sinf + 72 sin f ] 


i sin 1 


e"*^* sin 1 



14 



DMITRY CHELKAK AND STANISLAV SMIRNOV 





Figure 3. Local rearrangement at z\ the bijection between configura- 
tions. Without loss of generality, we may assume that the (reversed, i.e., 
going from 6"^) interface enters the rhombus z through the edge [^2^1] 
(the case [^114^2] is completely similar). There are two possibilities: ei- 
ther 7"^ (finally) leaves z through [uiiui] ("L", left turn) or through [^2^2] 
("R", right turn). In view of (12. 2p . the relative weights of configurations 
are \/2sini^, sin i^* ("L" pairs), and sin i^, ^f2^m\Q* ("R" pairs). 



where (y? denotes the total turn of the interface traced from to [n2U^i]. An easy 
trigonometric calculation using that 6 + 9* = 11/2 shows that 

[72 sin I + sin ^] - zsin^ = e*^ sin + e"'^^* [^2 sin | + sin ^] . 

Denoting the common value of the two sides by e^^ ■ F^{z) and observing that 1 ± z 
and e*^ ± e~*^*, we conclude that the first row ("L") describes the four projections 
of ■ F^{z) onto the lines M, e*^M, iM and e~*^*M, respectively. Multiplying by the 
common factor e~*2^ (which is always parallel to [i(wi — M2)]~2 ), we obtain the result 
for "L" pairs of configurations. The interchanging of 9 and 9* yields the result for "R" 
pairs. □ 

Remark 2.3. For Q G dVL^^ we define F^{Q so that /\2.5\) holds true (in this case, only 
two projections are meaningful, so F^{Q is easily and uniquely defined). Note that 
all interfaces passing through the half-rhombus ( intersect both its sides. Moreover, 
since the winding of the interface at ( is independent of the configuration chosen (and 
coincides with the winding of the corresponding boundary arc ) for topological reasons, 
we have 

F'iOWiriOrK Cednl, (2.6) 

where (see Fig. [HSj 

r(C) = «2(C)-^i(C), Ce(6V), «i,2(C)er, 

is the " discrete tangent vector" to dfli directed from to b^ on both boundary arcs. 
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Thus, we arrive at 



Discrete Riemann boundary value problem for (FK-case): The function 
is defined in ill and for each pair of neighbors Zo,Zi G fil, zq Zi, the discrete 



holomorphicity condition holds: 



Proj 



F\zo);[tiw-u)]-^ 



Proj 



F\z,);[iiw~u)]-^ 



(2.7) 



Moreover, F^ satisfies the boundary conditions Ii2.6\) and, since all interfaces pass 
through b^ , satisfies the normalization F^{b^) = lie F^{b^^) = (25)~2. 



2.2. Critical spin-Ising model. 



2.2.1. Definition of the model, holomorphic fermion, martingale property. Let f2^C<) 
be a simply-connected discrete domain composed of inner rhombi z G Int and bound- 
ary half-rhombi ( G dQ^^, with two marked boundary points a^,b^ G dfl^^, such that 
dfl^^ contains only "white" vertices (see Fig. HJA.) and there is no edge of F* breaking 
into two non-connected pieces. 
To each inner "white" vertex w G Int fif*, we assign the spin a{w) (+ or — ), thus 
obtaining the set of configurations (whose cardinality is 2*'^^"*'^r*)). We also impose 
Dobrushin boundary conditions assigning the — spins to vertices on the boundary 
arc (a^b^) and the + spins on the boundary arc (b^a^) (see Fig. HJA.). The partition 
function of the critical spin-lsing model is given by 

Z^ns.^s^^s) = [sin l9{b^)]-^ J2 Uw^^w2:a{w,)^a{w2)^^i^2, x^,^, = tan l9{z) , 

spin config. 

(2.8) 

where 9{z) is the half-angle of the rhombus U1W1U2W2 having center at z (i.e., tan 6yj^^^ = 
|w2 — Wi|/|'U2 — Mil). The first factor sin~^ doesn't depend on the configuration, and is 
introduced for technical reasons. 

Due to Dobrushin boundary conditions, for each configuration, there is an in- 
terface 7*^ running from to b^ and separating -|- spins from — spins. If F is not a 
trivalent graph, one needs to specify the algorithm of "extracting 7*^ from the picture" , 
if it can be done in different ways. Below we assume that, 

if there is a choice, the "interface" takes the left-most possible route (see Fig. 

With this choice the interface separates clusters of "+" spins connected through edges 
and clusters of "— " spins connected through vertices. Any other choice would do for 
the martingale property and eventual conformal invariance, as discussed in the Intro- 
duction. For example, one can toss a coin at each vertex to decide whether "+" or "— " 
spin clusters connect through it. Note that, drawing all the edges separating + spins 
from — spins, one can rewrite the partition function as a sum over all configurations 
zu of edges which consist of a single interface running from to b^ and a number of 
loops. Namely, 

^(n*;a*-^6*) = [sin 1^(6^)]"^ n«,i~«,2:[,«i;«,2] intersects 3;«,i«;2 • (2.9) 

ro={interface+loops} 
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Figure 4. (A) Ising model on isoradial graphs: discrete domain fi^ 
and a sample configuration (the partition function is given by fl2.8p ). 
By our choice of the "turning rule," loops and the interface 7*^ separate 
clusters of "+" spins connected through edges and clusters of "— " spins 
connected through vertices. To illustrate this, 7^ and loops are drawn 
slightly closer to "+" spins. The component of fl^ not "swallowed" by the 
path unshaded. The vertex ui is shaded since it is not con- 

nected to anymore. The vertices U2 and are shaded too since each 
of them is connected to the bulk by a single edge which contradicts our 
definition of connected discrete domains. Moreover, since the "interface" 
a'' 7^ -w 2 could arrive at each of these points only in a single way, our 
observable certainly satisfies boundary condition (12. lip at Z2 and z^, thus 
not distinguishing them from the other boundary points. (B) Two sam- 
ples of "interface pictures" composed from a number of loops and a single 
interface 7'' running from to z. To define the winding (7*^; w z) un- 
ambiguously, we draw 7"^ so that, if there is a choice, it turns to the left 
(for z G dfl^^, this corresponds to the edge-connectivity of "-I-" clusters). 



For z G Qa,, the holomorphic fermion is defined as (cf. (12. 3p . (12. 4p ^ 



7 117 windine(7 ; a -^2) 
^ l^J - ^in^;a^b^){^) - ~ ^-iwinding(a^^6^) ' ^^'^^^ 



As in (12. 9p . Z(^QS.a^^z«) denotes the partition function for the set of "interfaces pictures" 
containing (besides loops) one interface 7*^ running from to z (see Fig. |1]B), 
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for each configuration we count all weights corresponding to the drawn 
edges, including tan ^9 (a^) and [cos ^9{z)]~^ = [sin ^9{z)]^^ ■ta.n^9{z) 
for the first and the last. 

Then E^q^.^j^^^^) will stand for the expectation with respect to the corresponding 
probability measure. Equivalently, one can take the partition function, multiplying 
the weight of each configuration by the complex factor Q-^^^'^'^^'^sii -^z ) (^^\^[q[^^ fQ^- 
z G Intf2^, may be equal to one of the four different complex values a, ia, —a, —ia 
depending on the particular configuration). Finally, 

J-'^(b^) II (r(6'^))^^/^, where rih^) = W2{h^) —Wi(h^), is a normalizing 
factor that depends only on the structure of nearb^ (see Section Wl]) . 

Here and below, for C, G dVt^^, r{() = W2{C) — u!i{C) denotes "discrete tangent vector" 
to (9fi^ oriented counterclockwise (see Fig. IHA. for notation). For any ( G dQ^^, the 
winding(7''; -w () is fixed due to topological reasons, and so 

F'iOUriOrK (ednl\{a'}, while F\a') \\ i{r{a'))-'^ . (2.11) 

Remark 2.4 (Martingale property). For each z G Intfi^, -^(Q*\[a*7'5]-7'5 ^ 
martingale with respect to the growing interface {a^ = 7q, 7^ , 7^, ...) (till the stopping 
time when •y^ hits or z becomes separated from ¥ by the interface, see Fig. \^). 

Proof. It is sufficient to check that F^ has the martingale property when 7*^ makes one 
step. Let --yO,^ denote all possibilities for the first step. Then, 

^(Q*;a*^6*) = tan i6'(a'^) • ^(n«\[a'5a[];a[-*6'5) + ••• + ^(!^«\[a*a^]; afj^f,*) ' (2-12) 

and so 



[tan^9{a^)]^^Z(^^s.^s^f,s) [tan^9{a^)]-^Z(^n6.^s^f,s) 
Taking into account that the difference 

winding(a'^ b^) — winding(aL b^) = winding(a'' — )■ Ol), 

we obtain 



"""^ ^ln^;as,,i^) [tan|^(a^)]-i ■ Ee-t— ^(-^-^-^) 

and so on. On the other hand, counting "interface pictures" depending on the first 
step as in f l2.12p . we easily obtain 

[tanl9{a^)]~^Z(^ns.^s^,^ ■ Ee-^"'"'i'"s(^'''^'~^^) 
which gives the result. □ 
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Figure 5. Bijection between the sets {Q^; -w zq) and [Q^; Zi) 
(local notations zq i, 6*0,1, u and Wi are given in Fig. HP). In cases I-III, 
the winding (pi = winding (7^; w zi) is unambiguously defined by 
Lfo = winding(7Q; -w zq). In case IV, there are two possibilities: (fi is 
equal to either v^o - 27r + 6*0 + 61 (IVa, IVc) or v^o + 27r + 6*0 + 61 (IVb). 



2.2.2. Discrete boundary value problem for . 



Proposition 2.5. For each pair of neighbors zo,zi G fi^ separated by an edge {wiu), 
we have 



Proj 



F\zo);[iiw,-u)y 



Proj 



F\z,);[iiwi-u)]-^ 



(2.13) 



The proposition amounts to saying that F^ is spin holomorphic as in Definition 13.11 
below (see Fig. [Tp for notation). 

Proof. The proof is based on the bijection between the set of "interface pictures" 
{fl^; ~^ Zq) and the similar set {fl^; Zi), which is schematically drawn in Fig. |5l 

The relative contributions of the corresponding pairs to F^{zq) and F^{zi) (up to the 
same real factor) are given in the following table: 
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I 


[cosi^o]-'-e-t'^ 


[cos|^i]-i ■e-t(^+^"+''i) 


II 


[cos|^o]"^tani^i ■ e't^ 


[cosiei]-i-e-5{^"-+^o+ei) 


III 




[cos l^i]-! tan i^o ■ e-i(^-"+^"+^i) 


IV 


[cos|0o]"^tani^i ■ e^t'^ 


[cos|^^i]-itan|^^o-e-5(^±2.+eo+eo 



where ip = winding(7'^; -w zq) — winding(a^ -w 6^) + argr(6^). Note that 
e'i'^ II [wi-tyo]"^ e~i('^+^»+^i) II [w2-wi]-^, in cases I & II, 
e"^'^ II [wo-w^i]~i e-i('^+^»+^i) II [wi-W2]~K in cases III & IV. 
A simple trigonometric calculation then shows that the (relative) contributions to both 
projections FToi[F^{zj); — m)]"2] for j = 0, 1 are equal to 1, tani6'i, tan 16*0 and 
tan tan 16*1 in cases I-IV, respectively. □ 

Summing it up, we arrive at 

Discrete Riemann boundary value problem for (spin-case): The func- 
tion is defined in fi^ so that discrete holomorphicity /i2.13\) holds for every pair of 
neighbors zq, z\ G Furth ermore, F satisfies the boundary conditions Ii2.11\) and is 
normalized at . 

3. S-HOLOMORPHIC FUNCTIONS ON ISORADIAL GRAPHS 

3.1. Preliminaries. Discrete harmonic and discrete holomorphic functions 
on isoradial graphs. We start with basic definitions of the discrete complex analysis 
on isoradial graphs, more details can be found in Appendix and our paper |CS08j . 
where a "toolbox" of discrete versions of continuous results is provided. 

Let r be an isoradial graph, and H be defined on some vertices of F. We define 
its discrete Laplacian whenever possible by 

[A'H]iu) := ^^tan^^,- [if («,)-//(«)], (3.1) 

where /if(M) = XIm^^m sin 26'^ (see Fig. [T]B for notation). Function H is called 
(discrete) harmonic in some discrete domain fl^ if = at all interior vertices of 
f2f . It is worthwhile to point out that, on isoradial graphs, as in the continuous setup, 
harmonic functions satisfy some (uniform w.r.t. 6 and the structure of (}) variant of 
the Harnack's Lemma (see Proposition IA.4p . 

Let u\u] E; f2p) denote the harmonic measure of E G dQf^ viewed from u G 
Intf2p, i.e., the probability that the random walk generated by fl3.ll) (i.e. such that 
transition probabilities at u are proportional to tan^^'s) on F started from u exits fip 
through E. As usual, cu^ is a probability measure on 5f2p and a harmonic function of 
u. If Jlp is bounded, then we have 

H{u)= ^ u;'^(M;{a};^]^) ■if(a), wGlntl^^, 
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for any discrete harmonic in f2p function H. Below we use some uniform (w.r.t. 5 and 
the structure of V) estimates of from |CS08] . which are quoted in Appendix. 

Let H be defined on some part of F or F* or A = F U F* and z be the center of 
the rhombus Wif2f3f4. We set 



[d'H]{z) := \ 



Furthermore, let F be defined on some subset of <). We define its discrete 9^-derivative 
by setting 

[d'F]iu) = Y.^Ws+i-Ws)F{zs), M e A = F U F* (3.2) 



(see Fig. [Tp for notation when m G F). Function F is called (discrete) holomorphic in 
some discrete domain C <) if d^F = at all interior vertices. It is easy to check 
that A'^ = Ad^d^, and so d^H is holomorphic for any harmonic function H. Conversely, 
in simply connected domains, if F is holomorphic on <), then there exists a harmonic 
function H = F{z)d^z such that d^H = F. Its components H\r and H\r* are 
defined uniquely up to additive constants by 

H{V2) - H{Vi) = F {^{V2 + Vi)) ■ {V2 - Vi), V2 ~ Vi, 

where both fi,t>2 G F or both vi,V2 G F*, respectively. 

It is most important that discrete holomorphic (on <)) functions are Lipschitz 
continuous in an appropriate sense, see Corollary I A. 71 For the sake of the reader, we 
quote all other necessary results in Appendix. 

3.2. S-holomorphic functions and the propagation equation for spinors. In 

this section we investigate the notion of s-holomorphicity which appears naturally for 
holomorphic fermions in the Ising model (see (12. 7p . (I2.13P ). We discuss its connec- 
tions to spinors defined on the double-covering of <) edges (see [MerOlj ) and essentially 
equivalent to the introduction of disorder operators (see |KC71t IRC07j and the refer- 
ences therein). We don't refer to this discussion (except Definition 13.11 and elementary 
Lemma [3. 2 p in the rest of our paper. 

Definition 3.1. Let C be some discrete domain and F : fi^ — )■ C. We call 
function F strongly or spin holomorphic, or s-holomorphic for short, if for each 
pair of neighbors zq, z\ G ^2^, ~ the following projections of two values of F are 
equal: 



Proj 

or, equivalently. 



F{zo);[i{wi-u)]-^ =Proj F{zi) -[liwi-u)]-^ (3.3) 



F(^i)-F(zo) = -tiwi-u)5-' ■ iF{z^)-F{zo)), (3.4) 
where {wiu), m G F, i^i G F*, is the common edge of rhombi Zq, Zi (see Fig. \^). 
Recall that orthogonal projection of X on u satisfies 
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which we used above. 

It's easy to check that the property to be s-holomorphic is stronger than the usual 
discrete holomorphicity: 

Lemma 3.2. If F : — )■ C zs s-holomorphic, then F is holomorphic in Q^^, i.e., 
[d^F]{v) = for all v G Int Qi. 

Proof. Let v = u eT (the case v = w E T* is essentially the same) and Fg = F{zs) (see 
Fig. [T]B for notation). Then 

n n n 

-iY,i^s+i~Ws)F, = -tJ2{ws+i-u){F,-F,+i) = -5. ^(F,-F,+i) = 0. 

s=l s=l s=l 

Thus, [d^F]{u) =0. □ 

Conversely, in a simply-connected domain every discrete holomorphic function 
can be decomposed into the sum of two s-holomorphic functions (one multiplied by i): 

Lemma 3.3. Let be a simply connected discrete domain and F : C be 

a discrete holomorphic function. Then there are (unique up to an additive constant) 
s-holomorphic functions -Fi, -F2 : C such that F = Fi -\- iF2. 

Proof. Let denote the set of all oriented edges ^ = [^hC,w] of the rhombic lattice 
<0 connecting neighboring vertices G £ ^f*- For a function F : fi^ — )■ C, we 

define its differential on edges (more precisely, dF a, 1-form on the edges of the dual 
graph, but there is no difference) by 

dF{\uwi]) := F{zi) - F{zo) , dP-.n^^C, 

(see Fig. [T]C for notation). Then, a given antisymmetric function G defined anti- 
symmetrically on is a differential of some discrete holomorphic function F (uniquely 
defined on fi^ up to an additive constant) if for each (black or white) vertex u G flj^ 
(see Fig. [33 for notation when u eT) two identities hold: 

n n 

J2G{[uWs]) = and ^ G'([mw,])(w,-m) = 0. (3.5) 

s=l s=l 

Indeed, the first identity means that G is an exact form and so a differential of some 
function F, and the second ensures that F is holomorphic by (13.21) : 

n n n 

^G{[uws]){ws-u) = ^{F{Zs) - F{zs-i)){ws-u) = - ^F{Zs){Ws+l-Ws). 

s=l s=l s=l 

Note that identities (13.51) are invariant under the antilinear involution G 1— ?■ G", where 



G»(M) — G{[uw])-i{w-u)6~^ 

On the other hand, from (13.41) we see that F is s-holomorphic iff [dP]"^ = dP. Thus, 
the functions Pj defined by dPi = \{dF + {dPf) and dFg = j-idP - {dPf) are 
s-holomorphic and do the job. Uniqueness easily follows from (13. 3p : If = Fi + iF2, 
and both functions are s-holomorphic, then (13.31) implies that -^1(2:0) = F2{,zq), and 
uniqueness follows. □ 
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The spinor S{^) := [i{w^—u^)] 2 
is naturally defined on T ("con- 
tinuously" around rliombi and ver- 
tices). In particular, 







=: a, 


5(6) = 


-5(^6) 


= e'^'a, 


5(^3) = 


-5(6) 


= ia, 


SiU) = 


-5(6) 


■ iO* 

= le a 



Figure 6. Double covering T of T (= edges of ()■). 



Following Ch.Mercat [MerOlj . we denote by T the double- covering of the set T of 
edges <) which is connected around each z G <) and each v G A (see |Mer01] p. 209). A 
function S defined on T is called a spinor if it changes the sign between the sheets. The 
simplest example is the square root [i{w—u)Y'2 naturally defined on T (see Fig. E]). We 
say that a spinor S satisfies the propagation equation (see |Mer01j p. 210 for historical 
remarks and Fig. O for notation) if, when walking around the edges 6)---;6 of ^ 
doubly-covered rhombus, for any three consecutive edges spinor values satisfy 

^(e,+2) = {cose,)-^ ■ 5(6+1) - tan^^, ■ 5(6), (3.6) 

where 6j denotes the half-angle "between" 6 and 6+1 ; i-^-; % — if 3 odd, and 



I — 6*, if j is even, 
n 

is s-holomorphic then the real spinor 



Lemma 3.4. Let fi^ be a simply connected discrete domain. If a function F : — t- C 



Ff{[uw]) := Re F{zo) ■ [i{w-u)] 



Re 



F{zi) ■ [i{w-u)] 



(3.7) 



satisfies the propagation equation h3. 6)) . Conversely, if Ff is a real spinor satisfying 
hS. then there exists a unique s-holomorphic function F such that \3. 7|) holds. 



Proof. Note that (13. 6p is nothing but the relation between the projections of the same 
complex number onto three directions a = [z(ti'i— ui)]~2 , e**-5~^)a = [i{w2—ui)]~'2 and 
ia = [^(^^2 — ^2)]^^. Thus, s-holomorphicity of F implies (13. 6p . Conversely, starting 
with some real spinor Fy satisfying (13. 6p . one can construct a function F such that 

Proj[F(z) ; [i{w — u)]~^] = Ff{[uw]) ■ [i{w—u)]~^ 



(which is equivalent to (13. 7p ) for any z, and this function is s-holomorphic by the 
construction. □ 
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3.3. Integration of for s-holomorphic functions. 

Lemma 3.5. If F, F : — )• C are s-holomorphic, then [d^{FF)](v) E zM, v E Vl 



Proof. Let v = u E T (the case v = w E T* is essentially the same) and denote 
Fg := F{zs) (see Fig. [TB for notation). Using (13. 4p . we infer that d^F"^ satisfies: 



)Ft = - 



6 ■ 5^(F, + F,+i)(F,-F,+i) = 2iS ■ ImJ^FsFs+i E 



s=l 



□ 



Therefore, 4[d\FF)]{u) = [d^{F+FY]{u) - [d\F-F)^]{u) E iR. 

In the continuous setup, the condition Re[9G](2;) = allows one to define the 
function Im J G{z) dz (i.e., in simply connected domains, the integral doesn't depend 
on the path). It is easy to check that the same holds in the discrete setup. Namely, 
if is simply connected and G : — t- C is such that Re[(9''G'] = in fi^, then the 

discrete integral if = Im G{z)d^z is well-defined (i.e., doesn't depend on the path of 
integration) on both fip C F and Q^, C F* up to two (different for F and F*) additive 
constants. 

It turns out that ii G = F^ for some s-holomorphic F, then H = Im f {F(z)Yd^z 
can be defined simultaneously on F and F* (up to one additive constant) in the following 
way: 



H{u) - H{wi] 



26- 



Proj 



F{zj);[i{wi-u)] 2 



u 



Wi, 



(31 



where (uwi), -u G F, Wi G F*, is the common edge of two neighboring rhombi Zq,Zi E <) 
(see Fig. [IP), and taking j = 0, 1 gives the same value. 

Proposition 3.6. Let be a simply connected discrete domain. If F : fi^ C is 
s-holomorphic, then 

(i) function H : — )■ C well-defined (up to an additive constant) by 113. S\) : 
(a) for any neighboring f i, f2 G fif C F or Vi,V2 E fif, C F* the identity 

H{v2) - H{vi) = lm[{v2~vi){F{l{vi+V2))y] 

holds (and so H = Im j {F{z)fd^ z on both F and V* ); 

(Hi) H is (discrete) subharmonic on F and superharmonic on T* , i.e., 

[/S.^H]{u) ^ and [/S~^H]{w) ^ 
for all n G Int i7f C F and w E Int fif, C F*; 

Proof. (i),(ii) Let z be the center of the rhombus U1W1U2W2. For j = 1, 2 we have 
[H{U2) - Hiw,)] + [Hiw,) - Him)] 

= 2 |Re [[t{wj-U2)Y^^F{z)] f - 2 |Re [[i{wj-ui)f^F{z)] f 



'i{w,-U2)Y^^F{z) + [-tiw--u^)]y^F(^ 

'i{w,-Ui)]'/^Fiz) + [-i{w--W)]'^^'W) 
t{ui-U2){F{z)f - t{ui-U2)(F{zjf 



lm[{u2-ui){F{z)f] . 
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The computations for H{w2)~H{wi) are similar. 

(iii) Let m G F and set Fg := F{zs) (see Fig. [TB for notation). Denote 



Proj 



Fs][i{ws-u)] 2 



Proj 



Fs-i][i{ws-u)] 2 



Knowing two projections of Fg uniquely determines it's value: 



i{tse 



-ids 



ts+ie 



ids 



sin 9. 



and so 



Im 



tan 9s ■ 



-25 ■ Im 



-s=l 

" /j-2„-2t0. o+ + 1+2 „2ie. 



sin 6*0 



3iarg(?is-M) 



Let t= 



■ exp 



'-s=l 

i arg ( m)]^^ 



sin 6', 



, where G M and the argument of the square 



root changes "continuously" when we move around m, taking s = 1, ..,n. Then 



Im 



s=l 



™2^ i6s —J— ^ _[_ ™2 



sin 9. 



s=l * 

where "=f" is "— " for all terms except x„xi (these signs come from our choice of 
arguments). The non- negativity of the quadratic form Q*^"^ (for arbitrary 9i, ..,9n > 
with 9i-\-..-\-9n = tt) can be easily shown by induction. Indeed, the identity 



reduces the problem to the non- negativity of Q^^^'s. But, if a, /3, 7 > and = tt, 

then 



(3) 



I f n-1 , Cn 



^Xl, Xfi^l, Xy, 



sins ^ 



1 1 
sm2 a ■ sm2 7 



X — 



sm2 7 



1 

sin 2 a 



-\ 2 



■ z 



sin 2 a ■ sin 2 (3 sin 2 /3 ■ sin 2 7 

This finishes the proof for f = m G F, the case v = w G F* is similar. The opposite sign 
of [A^H]{w) comes from the invariance of definition (13. 3p under the (simultaneous) 
multiplication of F by z and the transposition of F and F*. □ 

Remark 3.7. As it was shown above, the quadratic form Qg"^. .^^(xi, .., x„) can be 

represented as a sum of n — 2 perfect squares Q^^\ Thus, its kernel is (real) two- 
dimensional. Clearly, this corresponds to the case when all the (complex) values F{zs) 
are equal to each other, since in this case [A^H] = by definition. Thus, 



(3.9) 



since both sides considered as real quadratic forms in xi,..,x„ (with coefficients of 
order 1 ) are nonnegative and have the same two-dimensional kernel. 
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3.4. Harnack Lemma for the integral of F^. As it was shown above, using 
for any s-holomorphic function F, one can define a function H = Im {F^{z)Yd^ z on 
both r and F*. Moreover, if|p is subharmonic while H\^^ is superharmonic. It turns 
out, that H, despite not being a harmonic function, a priori satisfies a version of the 
Harnack Lemma (cf. Proposition IA.4( i)). In Section [5.11 we also prove a version of 
the boundary Harnack principle which compares the values of H in the bulk with its 
normal derivative on a straight part of the boundary. 

We start by showing that and H\^, cannot differ too much. 

Proposition 3.8. Let w G IntQp. be an inner face surrounded by inner vertices 
Us G Int fip and faces Wg G fif*, s = l,..,n. If H is defined by liS. 8\) for some 
s-holomorphic function F : fl^ — )■ C, then 

max H{us) — H{w) ^ const ■ I H{w) — min H{ws 

s=l,..n \ s=l,..,n 

Remark 3.9. Since definition \3. 3\) is invariant under the (simultaneous) multiplica- 
tion of F by i and the transposition ofV and T* , one also has 

H{u) — min H(ws) ^ const ■( max H{us) — H{u)) 

s=l,..n s=l,..,n 

for any inner vertex u G Int VL^ surrounded by Wg G Int VLy* '^'^^ G fip, s = 1, .., n. 

Proof. By subtracting a constant, we may assume that m.ms=i^,,^nH{ws) = 0. Since 
H\-p, is superharmonic at w, it is non-negative there and moreover 

H{w) ^ const -H^Ws) for all s = 1, .., n. 

Since H\^{ws) ^ 0, by (13. Sp . it is non- negative at all points of F which are neighbors of 
Ws's. From subharmonicity of we similarly deduce that H{us) ^ const -if (ms+i). 
Therefore, 

H{us) X M := max H{us) for all s = 1, ..,n. 

s=l,..n 

We need to prove that K := M/H{w) ^ const. Assume the opposite, i.e. K ^ 1. 
Then, for any s = 1, ..,n, one has 



H(u.)-H{w) 

By (13.81) . these increments of H are derived from projections of F{zs) on two directions, 
and so 

|Proj[F(2:,); - ^ i^q(^j^-i\ 

|Proj[F(z,); [t{w ~ Us)]'^\ 
Here Zg is a center of the rhombus wUgWsUs+i, and we infer 



—1^ 



aig F{zs) mod TT = either a.Tg[{ws~w) 2] or a.Tg[{w—Ws) + 0{K 

Moreover, 6 ■ \F{zs)\'^ x H{us)—H{w) x M for all s = 1, .., n, if K is big enough. Since 
F{zs) and F{zs+i) have very different arguments, using (13. 9 p we prove the proposition: 

n 

M X 5 - ^ \F{zs+i)-F{zs)\^ X 6^ ■ \[A^H]{w) \ ^ const -if (w). □ 

s=l 
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Remark 3.10 (uniform comparability of and H^*). Suppose ^^|p, ^ on 
dQ^* o,nd hence, due to its superharmonicity, everywhere inside fl^, . Then, definition 
l[3. 8\) and Proposition \3.8\ give 



H^,{w) ^ H^{u) ^ const -H^^iw) 



for all neighboring u ^ w in Vt^ , where the constant is independent of and 6. 

Proposition 3.11 (Harnack Lemma for Im f^{F{z))^d^z). TakevQ G A = TUT* 
and let F : B^{uo, R) ^ C be an s-holomorphic function. Define 

1-5 

H ■.= lm / {F{z)yd^z -.ni^R 

by so that H ^ zn B{{uo, R). Then, 

H{vi) ^ const ■H{vq) for any Vi G B^{uo, |r). 

Proof. Due to Remark 13.101 we may assume that vo G F* while vi G F. Set M := 
max^g^a^^^ i^) _&(«). Since the function is subharmonic, one has 

M = H{vi) ^ H{v2) ^ H{v^) ^ . . . 

for some path of consecutive neighbors = {f i ~ f 2 ~ f 3 ~ • • • } C F running from 
Vi to the boundary 9-Bf(fo,-R). Take a nearby path of consecutive neighbors K^* = 
{wf ~ ~ It's ~ ...} C F* starting in (95^, (fo, |-R) and running to dB^*(vo, R). By 
Remark 13.101 H\^^ ^ const -M on K^,, and so 

H{vo) ^ const -M ■ uj\vo; /sT^.; 5^.(t;o, R) \ K^,) ^ const -M, 

since -f^|p, is superharmonic and the discrete harmonic measure of the path K^* viewed 
from vq is uniformly bounded from below due to standard random walk arguments 
(cf. [Ug08] Proposition 2.11). □ 

3.5. Regularity of s-holomorphic functions. 

Theorem 3.12. For a simply connected fi^ and an s-holomorphic F : — )■ C define 
H = {F{z)yd^z : fi^ C by (13. 8p in accordance with Proposition \3.(A Let point 
Zq G Int be a definite distance from the boundary: d = dist(zo; dQ^^) ^ const -6 and 
set M = max^^QS \H{v)\. Then 



Ml/2 

|F(^o)| ^ const (3.10) 



and, for any neighboring Zi ^ zq. 



MV2 



\F {zi)-F{zo)\^ const -^^-S. (3.11) 

Remark 3.13. Estimates 113. 10\} and i\3.11\} have exactly the same form as if H would 
be harmonic. Due to 113. 8\} and 113. 9\} . a posteriori this also means that the subharmonic 
function H\y and the superharmonic function H\y* should be uniformly close to each 
other inside fi^, namely H\ — H\ = 0{5M/d), and, moreover, \A^H\ = 0{5M/d^). 
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The proof consists of four steps: 

Step 1. Let By{zq; r) C F denote the discrete disc centered at zq of radius r. Then the 
discrete Li norm (as defined below) of the Laplacian of H satisfies 

ll^'^^IL;Bf,»:f., - ll^'^^l('')l''r(«) < oonst.M (3.12) 

and the same estimate for H restricted to T* holds. 

Proof of Step 1. Represent H on B^^Zq] d) as a sum of a harmonic function with the 
same boundary values and a subharmonic one: 

and 

i^harm = H, H^uh = OU dB^^ZQ] d) . 

Note that the negative function ifsub satisfies 

since the Green's function G{- ; u) in B^^zq; d) is negative. By the maximum principle, 
l-f^harml ^ M and SO l-ffsubl ^ 2M in By{zq-^ d). Therefore, 

const -Mrf^ ^ WH^nhWl, BUzo;d) ^ II ^''^sub || l ;B^(^o;|d) ' . \\G{- ■,u)\\^ dy 

Since A^H^nh = A^H, the inequality f l3.12p follows from the (uniform) estimate 

l|G'(- ;M)||i;Bf(^o;d) > const -ci^ for all u e B^{zo] fd) 
which we prove in Appendix (Lemma IA.8|) . □ 

Step 2. The estimate 

ll^ll2;B^(.o;i.) ^= I^WlVi(^) ^ const -Md (3.13) 

holds. 

Proof of Step 2. Since H = Im J^{F{z))'^d^z on both F and F*, it is sufficient to prove 

and a similar estimate for (9'^[i7|r.] . Represent H on B^^^zq; jd) as a sum: 
-f^lr = -f^harm + -f^sub, ^^-f^hami = 0, A^ifgub ^0 in -Br(^o; f c^) 

and 

^harm = H, H,^^ = OU dBf.{zo; fci). 
It follows from the discrete Harnack's Lemma (see Corollary \A.5\i and the estimate 

l-f^harml ^ M that 

|<9'^i^harm(^)| ^ COUSt -M/d 

for all z G Bf.{zo; jd) and hence 

ll^'^harm||i.5.(,^.l,) ^COnst-Mrf. 
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Furthermore, 



[d'H,^^]{z) = Yl [d'G]{z;u)[A'H]iu)4{u) 



ueInt_B^(zo;|d) 



where G{- ;u) ^ denotes the Green's function in B^{zo; ^d). We infer that 



4" 



The first factor is bounded by const -M (Step 1) and the second, as in the continuous 
setup, is bounded by const -d (see Lemma lA.9|) . which concludes the proof. □ 



Step 3. The uniform estimate 



Ml/2 

|F(2)| ^ const--^ (3.14) 



holds for all z G B^{zo; ^d). 

Proof of Step 3. Applying the Cauchy formula (see Lemma [A. 61 (i)), to the discs 

B^iz; 5Sk) C B^izo; irf), k : \d < 5Sk < {d, keZ , 
whose boundaries don't intersect each other and summing over all such k, we estimate 

d I r./ M ^ const \\^\\l;B'l{zo;U) 

f — . 

Thus, by the Cauchy-Schwarz inequality and f l3.13p 
const const ,, 

Ml/2 



Step 4- The estimate 



holds for all Zi ^ Zq. 



Ml/2 

|F(^i)-F(^o)| ^const-^-5. (3.15) 



Proof of Step 4- It follows from the discrete Cauchy formula (see Lemma lA.61 (ii)) 
applied to the disc B^{zq] ^d) and Step 3 for its boundary that there exist A,B E C 
such that, for any neighboring Zi ~ zq (see Fig. [T]C for notation), we have for both 
points the identity 

F(2;o,i) = Proj[A; Mo,i-'«J+Proj[5 ; Wo,2-u^iJ+0(e) = A+Proj[5-A; Wo^2-Wi\+0{e), 
where 

_ max^gaB^(^^.i^) \F{z) \ ■ 6 m^/^ ■ 6 

^~ d ^ rf3/2 ■ 

The definition of an s-holomorphic function stipulates that 

Proj[F(zo) - F{zi) ■ [i{w,-u)]-'2] = 

(for all Zi ~ zq), which implies B — A = 0(e), and hence F{zi) = A + 0(e) for 
2 = 0,1. □ 
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3.6. The "(r(2;))^" boundary condition and the "boundary modification 

trick". Throughout the paper, we often deal with s-holomorphic in functions 
satisfying the Riemann boundary condition 

m II (r(^))-^ (3.16) 

on, say, the "white" boundary arc Lf* C dfl^^ (see Fig. [7JA.), where r (2) := W2{z)—wi{z) 
is the "discrete tangent vector" to dfl^^ at Being s-holomorphic, these functions 
posses discrete primitives 

:= Im / {F\z)Yd^z. 

Due to the boundary condition fl3.16p . an additive constant can be fixed so that 

H^, = on L^,. 

The boundary condition for is more complicated. Fortunately, one can reformulate 
it exactly in the same way, using the following "boundary modification trick": 

For each half-rhombus UintWiW2 touching the boundary arc Lf , , we draw 
two new rhombi UintWiUiw and UintWU2'W2 so that the corresponding an- 
gles 61 = 02 are equal to \0 (see Fig. 

In general, these new edges (MintWi_2) constructed for neighboring inner vertices Mint, 
may intersect each other but it is not important for us (one can resolve the problem of 
a locally self-overlapping domain by placing it on a Riemann surface). 

Lemma 3.14. Let UintWiW2 be the half-rhombus touching L^* and ( = |(w2 + wi). 
Suppose that the function F^ is s-holomorphic in fi^ and F^{() \\ {w2 — Wi)^^. Then, 
if we set 

H^{U2) = H'^{ui) := H'^,{W2) = H'^,{wi), 
the function remains subharmonic at Mint ■ 

Proof. Definition (13. 8p says that 

i/^(Mint) - H^,{wr,2) = 25|Proj[F^(C) ; [^(«;i,2-«int)]-^]|' = 2(5cos2 f ■ \F\0\' 
because i{wi^2 — Uint) TT e^'^^{w2 — Wi) and F^{Q \\ {w2—Wi)~^. Therefore 

2 tan I ■ (if^(ui,2) - H^iuint)) = -25sin^ ■ |F^(C)|' 
and, if u denotes the fourth vertex of the rhombus UintWiuW2, 

tane ■ (H^iu) - H^iuint)) = tanO ■ Im[(F^(C))^(u-Mint)] = -25 sin ^ ■ {F^Ql"^- 

Thus, the standard definition of at u and the new definition at Mi^2 give the same 
contributions to the (unnormalized) discrete Laplacian at Mint (see (13.11) ). □ 

Remark 3.15. After this trick, an additive constant in the definition of can be 
chosen so that both 

H^, =0 on L^, and Hf. = on Z^, 
where L~ denotes the set of newly constructed "black" vertices near L^*. 
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(b'a% 



(bV) 





(B) 



(a'b') 




(A) 



Figure 7. "Boundary modification trick": (A) local modification of 
the boundary; (B) FK-Ising model: example of a modified domain; 
(C) spin-Ising model: example of a modified domain. 



4. Uniform convergence for the holomorphic observable in the 

fk-ising model 



Recall that the discrete holomorphic fermion F^{z) 



F^QS.^s^hS){z) constructed in 



Section 12.11 satisfies the following discrete boundary value problem: 

(A) Holomorphicity: F^{z) is s-holomorphic in fi^. 

(B) Boundary conditions: F^(() \\ {t(Q)~^ for ( E dVl^, where 

t{0 = MO-MO, Ce{a'b'), wi,2(C)er*, 
r(C) = «2(C)-«i(C), Ce{b'a'), Mi,2(C)er, 

is the "discrete tangent vector" to dfl^^ directed from to on both 
arcs (see Fig. \^). 



(4.1) 



(C) Normalization at b^: F^{b^) = ReF^(6^ 



(25)- 



Remark 4.1. For each discrete domain {Q^^; ,b^) , the discrete boundary value prob- 
lem (A)&(B)&(C) has a unique solution. 

Proof. Existence of the solution is given by the explicit construction of the holomorphic 
fermion in the FK-Ising model. Concerning uniqueness, let F^ and F^ denote two 
different solutions. Then the difference F^ := -F/ — F2 is s-holomorphic in thus 
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z is well-defined (see Section 13. 3[ especially (13. 8p ). Due to 
condition (B), is constant on both boundary arcs (a^b^) C T* and (b^a^) C T. 
Moreover, in view of the same normalization of near one can fix an additive 
constant so that H^. = on (a^b^) and if^ = on (b^a^). 

The "boundary modification trick" described in Sect. 13.61 provides us the slight 
modification of fi^ (see Fig. 03) such that the Dirichlet boundary conditions = 0, 
ifp* = hold true everywhere on dQ\. Using sub-/super-harmonicity of on T/T* 
and dSSD, we arrive at ^ if^ ^ ^ in . Thus, = and Ff = F|. □ 

Let ^ 

H' = ■■= / iFlnS;aM>s)i^)Yd'z. (4.2) 

It follows from the boundary conditions (B) that is constant on both boundary arcs 
(a^b^) C r* and (b^a^) C F. In view of the chosen normalization (C), we have 

Remark 4.2. Due to the "boundary modification trick" (Section \3.6\} . one can fix an 

additive constant so that 

Hl^ = on {a%^)f, H^^, = on (a^b^), 

H^ = l on (b^a^), H^, = 1 on {b^a%,, ^ ' ' 

where {a^b^)j^ (and, in the same way, {b^a^)j^^) denotes the set of newly constructed 
"black" vertices near the "white" boundary arc (a^b^) (see Fig. \7p). 

Let f^{z) = f^Qs.g^s f,s-^{z) denote the solution of the corresponding continuous 
boundary value problem inside the polygonal domain Q^: 

(a) holomorphicity: is holomorphic in ; 

(b) boundary conditions: /''(C) || (t(C))^^ for C ^ dVl^ , where t{Q 
denotes the tangent vector to d^l^ oriented from to b^ ( on both arcs ); 

(c) normalization: the function = h^s i^s '■= Im Jif^iOYd^ is 
uniformly bounded in fl^ and 

Note that (a) and (b) guarantee that is harmonic in and constant on both 
boundary arcs {a^b^), (b^a^). In other words, 

f^ = V2tdh^, = u{-;b^a^;Q^), 

where u denotes the (continuous) harmonic measure in the (polygonal) domain Q^. 
Note that dh^ 7^ in since is the imaginary part of the conformal mapping 
from onto the infinite strip (—00, 00) x (0, 1) sending and b^ to =Foo, respectively. 
Thus, f^ is well-defined (up to the sign). 

Theorem 4.3 (convergence of FK-observable). The solutions F^ of the discrete 
Riemann-Hilbert boundary value problems (A)&(B)&(C) are uniformly close in the bulk 
to their continuous counterpart f^ defined by (a)&(b)&(c). Namely, for all < r < R 
there exists e{6) = e{S, r, R) such that for all discrete domains (f^^; a^, b^) and z^ G fi^ 
the following holds true: 
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ifB{z\r) cn^ C B{z\R), then \F\z^) - f{z^)\ ^ e{5) ^ as 5 ^ 

(for a proper choice of f^'s sign), uniformly with respect to the shape of and . 

Remark 4.4. Moreover, the sign of f^ is the same for, at least, all 'z^ lying in the 
same connected component of the r -interior of . 

Proof. Assume that neither f^ nor — /"^ approximates well, and so for both signs 
\F^{z^)±f^{z^)\ ^ £o>0 for some sequence of domains fl^, 5^0. Applying translations 
one can without loss of generality assume z^ = for all 5's. The set of all simply- 
connected domains Q : B{0,r) C C -8(0, R) is compact in the Caratheodory topology 
(of convergence of conformal maps germs). Thus, passing to a subsequence, we may 
assume that 

[Q^; a\ h^) ^ (^]; a, 6) as 5 ^ 

(with respect to = z^). Let h = h(^Q.a,b) '■= ba; Q). Note that ^ h as S 0, 
uniformly on compact subsets of Q, since the harmonic measure is Caratheodory stable. 
Moreover, 

if^)^ = 2idh^ ^ f^ = 2idh as 5 0. 
We are going to prove that, at the same time, 

and {F^f ^ f^ as 5 ^ 0, 

uniformly on compact subsets of Q, which gives a contradiction. 

It easily follows from (14.31) and the sub-/super-harmonicity of on F/F* that 

^ ^ 1 everywhere in fi^. 

In view of Theorem 13.121 this (trivial) uniform bound implies the uniform bounded- 
ness and the equicontinuity of functions F^ on compact subsets K of Q. Thus, both 
{H^} and {F^} are normal families on each compact subset of Q. Therefore, taking a 
subsequence, we may assume that 

F^ ^F and ^ H for some F : n C, H : n ^ R, 

uniformly on all compact subsets of fl. The simple passage to the limit in (14.21) gives 

H{V2) - H{v,) = lm [ {FiOfdC 

J[vi;V2] 

for each segment [vi; V2] C Q. Thus, F"^ = 2idH. Being a limit of discrete subharmonic 
functions H^, as well as discrete superharmonic functions H^*, the function H should 
be harmonic. The sub-/super-harmonicity of on F/F* gives 

00^ { ■ ; (6V)p,; fi^.) ^ H^, ^H^^uj^i-; (6V)r; ^r) in ^A, 

where the middle inequality holds for any pair of neighbors w G F* , u G F due to (13. 8p . 
It is known (see |CS08 J Theorem 3.12) that both discrete harmonic measures uj^{-) (as 
on F*, as on F) are uniformly close in the bulk to the continuous harmonic measure 
uj{-) = h. Thus, =^ h uniformly on compact subsets of f2, and so F"^ = 2idh. □ 
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Proof of Remark \4.4\ Consider simply-connected domains {Q^; with 
lying in the same connected components of the r-interiors Qf. Assume that we have 
\F\z^) - f^{z^) \ but \F^(z^) + f^(z^) \ ^ as 5 ^ 0. Applying translations to 
and taking a subsequence, we may assume that 

in^;a^,b^;z^)^{n;a,b;I) w.r.t. = /, 
for some F connected with inside the r-interior Qr oi Q. As it was shown above, 



F^ ^ F and =^ / uniformly on 



where either F = f or F = —f (everywhere in Q), which gives a contradiction. □ 

5. Uniform convergence for the holomorphic observable in the 

spin-ising model 

For the spin-Ising model, the discrete holomorphic fermion F^{z) = F^^s.^^s i,s-^{z) 
constructed in Section 12.21 satisfies the following discrete boundary value problem: 
(A°) Holomorphicity: F^{z) is s-holomorphic inside fi^. 
(B°) Boundary conditions: F^(Q \\ (r(^))^^ for all ( E except 
at , where t{Q = W2(C)— 'U^i(C) ^■^ the "discrete tangent vector" to dQ^^ 
oriented in the counterclockwise direction (see Fig.\^). 
(C°) Normalization at b^: F^{b^) = J^^{b^), where the normalizing 
constants J^^{b^) \\ (r(6''))^2 are defined in Section [3T71 

Remark 5.1. For each discrete domain (fi^; a*^, 6*^), the discrete boundary value prob- 
lem (A° )&(B° j&fO' ) has a unique solution. 

Proof. Existence is given by the holomorphic fermion in the spin-Ising model. Con- 
cerning uniqueness, let F^ denote some solution. Theni7^ = J {F^{z)fd^ z is constant 
on (9fif*, so either F^{a^) \\ (r(a'^))-i/2 F^{a^) \\ (-r(a'^))~i/^ In the former case, 
using the "boundary modification trick" (Section 13.61) . we arrive at = on both 
and dQ~. Then, sub-/super-harmonicity of on T/T* and (13. 8p imply that 
^ ifp ^ ifp* ^ in Q^. Therefore, = and F^ = 0, which is impossible. Thus, 
F^{a^) II (— r(a'^))~^/^ (for the holomorphic fermion this follows from the definition). 
Let Fi^2 be two different solutions. Denote 

F'{z) := (-r(a^))^ ■ [F/(a^)F|(z) - F^{a')Ft{z)]. 

Then, F^ is s-holomorphic and, since (—T(a^))i Ff 2(0,^) £ satisfies the boundary 
condition (B°) on dVl^ \ {a^}- Moreover, we also have F^{a^) = || (r(a^))~^. Arguing 
as above, we obtain F^ = 0. The identity Ff = F2 then follows from (C°). □ 



Let 



Remark 5.2. Using Section one can fix an additive constant so that 

Hft = everywhere on dfl^. 



Hf = everywhere on d^l~ except a 



"out' 



(5.1) 
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where dVL~, denotes the modified boundary ( everywhere except ) and a^^^ is the original 
outward "black" vertex near (see Fig. \7jp). Then, ^ H^, ^ everywhere in . 

5.1. Boundary Harnack principle and solution to (A°)&(B°) in the discrete 
half-plane. We start with a version of the Harnack Lemma (Proposition 13. lip which 

compares the values of = Im y{F^ f{z)d^ z in the bulk with its normal derivative 
at the boundary. 

Let R{s,t) := (— s; s) x (0; t) C C be an open rectangle, R^{s,t) C T denote its 
discretization, and L^{s), U^{s,t) and V^{s,t) be the lower, upper and vertical parts 
of the boundary (see Fig. EJA.). 

Proposition 5.3. Let t ^ 6, be an s-holomorphic function in a discrete rectangle 
i?^(2t,2t) satisfying the boundary condition (B°) on the lower boundary L^^{2t) and 

= Im J^{F^{z))'^d^z be defined by i\3. 8\} so that H = on Lf . and H ^ everywhere 
in i?p*(2i, 2i). Let E be the boundary vertex closest to 0, and G T* denote the 
inner face (dual vertex) containing the point c = it. Then, uniformly in t and 6, 

Proof. Let t ^ const -5 (the opposite case is trivial). It follows from Remark 13.101 
and Proposition 13.111 that all the values of on U^(t, |t) are uniformly comparable 
with H(c^). Then, the superharmonicity of -f^'^jp* and simple estimates of the discrete 
harmonic measure in i?f.(t, |t) (see Lemma [Ojl give 

H^{bi.) ^ io^ibf^.; f/^,; i?^*(t, ■ min^^g^^^^ H^{w^) ^ const -d/t ■ H\c^), 

where 6f, G F* denotes the inner dual vertex closest to (see Fig. [HJA.). Therefore, 
^ const ■iJ'^(c'^) /if: for a neighbor of b^. Due to the s-holomorphicity of F^, this is 
sufficient to conclude that 

|F^(6^)P ^ const ■iJ^(c*)/t. 

On the other hand, since = on L^{2t), one has ^ const -H^i^c^) everywhere 
in R^ity-^t) (the proof mimics the proof of Proposition 13.111 if H^{v) S> H^{c^) at 
some V G -Rf (t, \t), then, since = on L^{2t), the same holds true along some path 
runnmg from V to U^{2t, 2t)UV\2t, 2t), which gives a contradiction). Thus, estimating 
the discrete harmonic measure in Rr{t,\t) from the inner vertex 6f G F closest to b 
(see Fig. [HJA.), we arrive at 

H\bi) ^ J{b^^- U V^- Ri{t, it)) ■ max„ig^*ui?* H\u^) ^ const -d/t ■ H\c^). 

Since i/^(6f) x 5 ■ \F^{b^)\^, this means \F\b^)\'^ < const ■H^{c^)/t. □ 

Now we are able to construct a special solution J-'^ to the discrete boundary value 
problem (A°)&(B°) in the discrete half-plane. The value J^^{b^) will be used later on 
for the normalization of the spin-observable at the target point b^ . 



Theorem 5.4. Let M.^ denote some discretization of the upper half-plane (see Fig. ) 

1^ 



Then, there exist a unique s-holomorphic function : Hi — )• C satisfying boundary 
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conditions J^^{C) II (t(C)) ^ for ( G 9EI^, such that 

J^\z) = 1 + 0(5^ . (Imz)-5), 
uniformly with respect to <)'^. Moreover, x 1 on the boundary 9EI^. 

Remark 5.5. // dU^ is a straight line (e.g., for the proper oriented square or trian- 
gular/hexagonal grids), then = 1 easily solves the problem. 

Proof. Uniqueness. Let Ff, F| be two different solutions. Clearly, := F^^ — F| is 
s-holomorphic and satisfies the same boundary conditions on 9EI^. Thus we can set 

:= Im j\F^{z)fd^z, where = on both fflf* and dM^ (see Section ESD- 
Since {F^{z))'^ = 0{6 ■ (lmz)~^), the integration over "vertical" paths gives 

H^{v) = 0{6 ■\og{6~^lmv)) as Imt> — )■ oo, 

so grows sublinearly as Imt; — )■ oo which is impossible. Indeed, using simple esti- 
mates of the discrete harmonic measure (see Lemma [A. 30 in big rectangles R{2n,n), 
n — )■ cxD, and sub-/super-harmonicity of on F/F* together with the Dirichlet bound- 
ary conditions on the boundary dMj^, we conclude that 

Hf.{u) ^ lim 0{6 ■ log((5^^n)) • {lmu + 26)n~^ = for any u G H^, 

n—^oo 

and, similarly, H^,{w) ^ for any w G Hp,. Thus, = 0, and so F^ = F2. 

Existence. We construct J-'^ as a (subsequential) limit of holomorphic fermions in 
increasing discrete rectangles. Let 6 be fixed, E (} denote the closest to boundary 
vertex, and i?^ denote discretizations (see Fig. [H]A-) of the rectangles 

Rn = i?(4n,2n) := {-An; An) x (0;2ra). 

Let F,^ : i?^ ^ — t- C be the discrete s-holomorphic fermion solving the boundary value 
problem (A°)&(B°) in with being the discrete approximations of the points 2m. 
For the time being, we normalize F^ by the condition 

\F'{bi)\ = i. 

Having this normalization, it follows from the discrete Harnack principle (Proposi- 
tions [3TTT] and [5]3]) that X n everywhere near the segment [—2n + in;2n + in]. 
Moreover, since if = on the lower boundary, one also has H ^ const -n everywhere 
in the smaller rectangle i?^(2n, n) (the proof mimics the proof of Proposition I3.11j ). 
Thus, estimating the discrete harmonic measure of U^{2n, n)r U y/(2n, n) in i?f (2n, n) 
from any fixed vertex G and using the subharmonicity of H\^, one obtains 

Im + 26 

H^iv^) ^ const -n ^ const ■{lmv^ + 25), 

if n = n{v^) is big enough. Moreover, since -f^'^jp, is superharmonic, one also has the 
inverse estimate for near the imaginary axis i]R+: 

H^{v^) ^ const ■{lmv^ + 26), if | Ret;^| ^ 6. 

Further, Theorem 13.121 applied in (Re 2; — | Im 2;; Re z + | Im 2;) x (| Im z; | Im z) gives 
l-^n(^'')l ^ const for any z^ G Hi, if n = n{z^) is big enough. 
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Note that there are only countably many points G and G H^. Since 
for any fixed vertex the values H^{v^) and F^{z^) are bounded, we may choose a 
subsequence n = — )• oo so that 

Hliv^) H\v^) and F^(z^) ^ J^\z^) for each G and G H^, 

It's clear that J"-^ : ^ C is s-holomorphic, = Im /^(J'''(z))2c/'^z ^ 0, and J""^ 
satisfy the same boundary conditions as H^, F^, and J-'^{h^) = 1. Moreover, 

n\v^) = 0{lmv^ + 25), T\z^) = 0(1) uniformly in B.\ 

and n^v^) X {lmv^ + 26) for t;'' near «+. ^^'^^ 

Now we are going to improve this estimate and show that, uniformly in H^, 
H^{v) = /i ■ (Imf + 0(5)) for some fi > 0. For this purpose, we re-scale our lattice and 
functions by a small factor e — )■ 0. Let 

v'^ := ev^ , ^""^ := ez^ , and T^^^ft;""^) := eUHv^) , J^^'^f^"^) := J^^(/). 



Note that the uniform estimates (15 ■2p remains valid for the re-scaled functions. There- 
fore, Theorem 13.121 guarantees that the functions and J-'^'' are uniformly bounded 
and equicontinuous on compact subsets of H, so, taking a subsequence, we may assume 

'H^^{v) =^ h{v) uniformly on compact subsets of H. 

Being a limit of discrete subharmonic functions H^^ as well as discrete superharmonic 
functions Hf^i, the function h is harmonic. Moreover, it is nonnegative and, due to 
(15. 2p . has zero boundary values everywhere on M. Thus, h{v) = fiv for some /i > 
(the case /i = is excluded by the uniform double-sided estimate of Ti'^^ near iM+). 
Thus, for any fixed s ^ t > one has 

V.^^ ^ /it on [—s+it; s+it] as e — 0. 

For the original function H^, this means 

1-L^(v^) = (/i + Ofe_>oo(l)) ■ Imf'' uniformly on U^{ks, kt) as k = e^^ — > oo. 

Estimating the discrete harmonic measure in (big) rectangles R^{ks,kt) from a 
fixed vertex (see Lemma [A. 3^ and using subharmonicity of i7 |p and superharmonic- 
ity of -f^lp,, we obtain 



kt \ {ks 



If I ■ kt 

(/i + Ofc^oo(l)) ■kt + 0[ ] ■ 0{kt), 



where the 0-bounds are uniform in v^, t and s, if A; = k{v^) is big enough. Passing to 
the limit as — >■ oo, we arrive at 

n\v^) = ■ (Im/ + 0{6)) + 0{\v^\ ■ t^s^), 

where the 0-bound is uniform in and t, s. Therefore, 

n\v^) = ^ ■ {\mv^ + 0(5)) uniformly in . 
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(A) (B) 

Figure 8. (A) Discretizations of the upper half-plane EI and the rec- 
tangle R{s, t) = (— s; s) X (0; t). Boundary points a^, 6'' G <) approximate 
the points it and 0, while G F* approximate the center ^it. We de- 
note by L^{s), U^{s,t) and V^{s,t) the lower, upper and vertical parts 
of dR^{s, t), respectively. (B) To perform the passage to the limit under 
the normalization condition at 6 , we assume that D Rl{s,t) and 
dfl^^ \ {a^} D L^^{s). For 6 small enough, the discrete harmonic measure 
from of any path = K^^ going from to is uniformly bounded 
from below. We can similarly use any point lying on the "straight" 
part of the boundary for the (different) normalization of the observable. 



It follows from (15.21) and Theorem 13.121 that both J^^ and (J^^)'^ are uniformly 
Lipschitz in each strip (3 ^ lm( ^ 2/3 with the Lipschitz constant bounded by 0(/3~^). 
Taking some v G near z and v' G such that \v' — f | x 5^^'^{lmzY^'^, we obtain 

H\v') - n\v) = Im / (^-^(0)'^ = lm[{J^\z)f{v'-v)] + O 

J[v,v'] 

i.e., lm[{T^{z)f{v'-v)] = /i ■ lm{v'-v) + 0{5) for all v'. Thus, 

(jr<5(^))2 = jj + 0{6^ ■ (Imz)-^) uniformly in H^. 

Since J^^ is Lipschitz with the Lipschitz constant bounded by 0((lmz)~^) (see above), 
this allows us to conclude that 

±J^^{z) = fi^ + 0{6^ ■ {lmz)~^) uniformly in 

(for some choice of the sign). Thus, the function := i/i'^/^J-""^ satisfies the declared 
asymptotics and boundary conditions. Moreover, \J'^(b^)\ = yU^^/^ x L Since such a 
function J^^ is unique, all other values \J^^\ on the boundary dUi are x 1 too. □ 
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5.2. Main convergence theorem. To handle the normahzation at of our discrete 
observable, from now on we assume that, for some s,t > 0, 

contains the discrete rectangle R^{s,t), 
dQ^^ \ {a^} contains the lower side L^^{s) of (5.3) 
and is the closest to vertex of dfl^^ {see Fig.^. 

Some assumption of a kind is certainly necessary: one can imagine continuous domain 
with such an irregular approach to b, that any approximation is forced to have many 
"bottlenecks," ruining the estimates. 

Let f^{z) = f^^s.^s f,s-j{z) denote the solution of the following boundary value 

problem inside the polygonal domain fl^ (here the tilde means that we slightly modify 
the original polygonal domain fl^ near b^ , replacing the polyline L^{s) by the straight 
real segment [—s; s], cf. |CS08j Theorem 3.20): 

(a°) holomorphicity: f^ is holomorphic in fl^ ; 

(b°) boundary conditions: f^{() \\ (t(C))^^ for ( G dQ^ , where t{() 
is the tangent to dQ^ vector oriented in the counterclockwise direction, 
f^ is bounded away from ; 

(c°) normalization: the function = hf^g i,s '■= Ini /(/'^(C))^'^C 
nonnegative in Vl^ , bounded away from , and 

As in Section m (a°) and (b°) guarantee that is harmonic in fl^ and constant 
on dfl^. Thus, 

f' = V2tdJ?, where h' = P(n.;,.,o), 

is the Poisson kernel in having mass at and normalized at 0. In other words, is 
the imaginary part of the conformal mapping (normalized at 0) from onto the upper 
half-plane H sending and to oo and 0, respectively. Note that dh^ ^ everywhere 
in Q^, thus f^ is well-defined in up to a sign, which is fixed by f^{0) = +1. 

Theorem 5.6 (convergence of the spin-observable). The discrete solutions of the 
Riemann-Hilbert boundary value problems (A° )&(C^ ) are uniformly close in the 
bulk to their continuous counterparts f^, defined by (a° )&(b° )&(c° ). Namely, there 
exists e{5) = e{6,r, R, s,t) such that for all discrete domains {fl^;a^,b^) C B{0,R) 
satisfying Ii5.3\) and for all e fl^^ lying in the same connected component of the 
r -interior of Vt^ as the neighborhood of b^ (see Fig. \^ the following holds true: 

\F\z^) - f\z^)\ e{6) ^ as 6-^0 

(uniformly with respect to the shape of Vt^ and the structure of ()^). 

Moreover, it is easy to conclude from this theorem that the convergence also 
should hold true at any boundary point such that has a "straight" part near . 
Namely, as in (15. 3p . let (see Fig. |8j3) 

G dVL^^ he a boundary point, the boundary dVL^^ is "straight" near d^ 
and oriented in the (macroscopic) direction Td : \Td\ = 1, i.e.. 
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fi^ contains the discretization R^{s,t) of the rectangle d + Td- R{s,t) and 
dil^^\{a^} contains the discretization L^^(s) of the segment d-\-rd-[—'s;J\. 

Further, let J^^ denotes the solution of the boundary value problem (A°)&(B°) in 
the discrete half-plane {d + ■ H)^ which is asymptotically equal to (r^)^^^^ (again, 
= (r^i)^^/^, if one deals with, e.g., the properly oriented square grid). 

Corollary 5.7 (convergence of spin-observable on the boundary). If centers of 
R^^{s,t) and R^^(s,t) are connected in the r-interior of C B{0,R), then 

\F\d^)-[{rdy/^P{d^)]-f^{d)\^e{6)^0 as 6-^0 

(one should replace both Ly*{s) and i^f. (s) by the corresponding straight segments to 
define properly the value of the "continuous" solution f^{d) on the boundary). 

Proof of Corollary 15.71 Let be a discrete approximation of the point c := and 

be a discrete approximation of the point c := d+ \iTdt. It follows from the discrete 
Harnack principle (Proposition 15.31 and Proposition IA.4p that 

\F\d^)\ X \H\^)\ X \H\c^)\ X \F\h^)\ X 1, 

uniformly in VL^ and 5, if all the parameters r, i?, s, t, s, t are fixed. Denote by the 
discrete observable F^ renormalized at d^: 

■ F^(rf^) ' 

and by f^ the corresponding continuous function renormalized at d: 

75 _ (''rf) ^^"^ rS 

^ f^d) ■ ^ • 

Again, \f\d)\ x \h\Z)\ x \h\c)\ x |/'5(0)| = 1 due to the Harnack principle. More- 
over, since is equal to the imaginary part of the conformal mapping from VL^ onto 
the upper half-plane, the Koebe Distortion Theorem gives 

\f\c^)\^ = 2\dh\c^)\^\h\c^)\^l. 

One needs to prove that the ratio 

F^{d^) _F^ p _ F\c^) f{c^) 

[{Tdy/^^^{d^)] ■ f\d) ~¥^'f~ f\c^) ' F\(^) 

is uniformly close to 1. This follows from Theorem 15. 6 [ since F^[c^) is uniformly close 
to f\c^), F\c^) is uniformly close to P{c^), and \f{c^)\ x \f\c^)\ x 1. □ 

Proof of Theorem 15.61 Assume that 

|fV)-/'(^')I^^o>o 

for some sequence of domains il^ with 5 — )■ 0. Passing to a subsequence, we may 
assume that — )■ z. The set of all simply-connected domains B{z,r) C f2 C B{0,R) 
is compact in the Caratheodory topology, so, passing to a subsequence once more, we 
may assume that 

{Q^; a^, b^) {Q; a, b) with respect to / ^ G i7 as 5 -> 0. 
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Note that Q D R{s, t) = {-s; s) x (0; t), dQ D [-s; s], and 6^ -> 6 = 0. Let h = h(yi.a,b) 
be the continuous Poisson kernel in Vl having mass at a and normahzed at (i.e., the 
imaginary part of the properly normalized conformal mapping from onto H). Then, 

/i'' =^ /i as 5 0, 

uniformly on compact subsets of f2, since this kernel can be easily constructed as a 
puUback of the Poisson kernel in the unit disc. Moreover, it gives 

f = V2idh^ ^ f = V2idh as 5 ^ 

uniformly on compact subsets of Q (here and below the sign of the square root is chosen 
so that f^{0) = /(O) = +1). We are going to prove that, at the same time, 

and ^ V2idh as 5 ^ 0. 

We start with the proof of the uniform boundedness of away from . Denote 
by c := ^it the center of the rectangle R{s,t) and by G P* the dual vertex closest 
to c (see Fig. 03)- Let d>0 he small enough and 7^ C B{ad, ^d) be some crosscut in 
fl separating a from c in fl. Further, let L^^ G fl^ H dB{ad, Sd) be an arc separating 
from in (such an arc exists, if S is small enough), and denote the connected 
components of fl^ \ L^^, containing c^. 

The Harnack principle (Propositions 15.31 and 13. lip immediately give 

H^{c^) X 1 uniformly in 6, 

if s and t are fixed. Let 

M|,:=max{i7^(w^), e i^lM- 

Because of the subharmonicity oi H\^, = Hf.{uf^) ^ Hf.{ul) ^ H^{u^^) ^ . . . for 
some path of consecutive neighbors = {uq ~ ~ ~ . . . } C P. Since the 
function vanishes everywhere on dQ~ except a^, this path necessarily ends at a^. 

Taking on the dual graph a close path Kf^, = {wq ~ iff ~ ^2 ~ . . . } C P* starting 
near Mq and ending near a^, we deduce from Remark 13. 101 that 

H^.iwi)^ const -M^,. 

Then, 

H^4c^) ^ J{c^; K^.;ni. \ ^r*) " const -M^^ ^ const((n; a), d) ■ M^^, 

since H\^^ is superharmonic and 

u\c^] K^.;^^, \Kl,)^ \u{c^- Kl,-Vt^ \K^t')> const((fi; a), d) > 

for all sufficiently small 5's (see Fig. HP and |CS08j Lemma 3.14). 

Thus, the functions are uniformly bounded away from . Due to Theo- 
rem [n]T2l we have 

= 0(1) uniformly on compact subsets of ^2. (5.4) 

Moreover, using uniform estimates of the discrete harmonic measure in rectangles 
(Lemma IA.3P exactly in the same way as in the proof of Theorem 15. 4[ we arrive at 

H\v) = 0{lmv), F^ = 0(1) uniformly in R^^s, ^t). (5.5) 
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Taking a subsequence, we may assume that 

and ^ H for some F : ^ C, H : n ^ R, 

uniformly on all compact subsets of Q. The simple passage to the limit in ( 14. 2 p gives 
H{v2)-H{vi) = Im /j^^.^^j(F(C))2dC, for each segment [vi;v2] C Q. So, F^ = 2idH, and 
it is sufficient to show that H = P(Q-a,b) ■ Being a limit of discrete subharmonic functions 
ifp, as well as discrete superharmonic functions -ffp*, the function H is harmonic. The 
next step is the identification of the boundary values of H. 

Let M e r2 and c? > be so small that u G ilf^. Recall that the functions H^\^ 
are subharmonic, uniformly bounded away from a^, and = on the (modified) 
boundary d^l~, except at . Thus, the weak Beurling-type estimate of the discrete 
harmonic measure (Lemma IA.2p easily gives 

H{u) = ^mH^iu) ^ const ■ im ... . Lsf o.^^ 

5->o (5-s>o L distf^a [u ] oBf,[ad, 3d)) 

^ const(r2, d) ■ (dist(M ; dVl-^d \ B^ad, 3d))'^ for all u G fis^ 

(since, if 6 is small enough, u G ^Id)- Thus, for each d > 0, H{u) — )■ as m — )■ 9f2 inside 
VL^d) i-e., = on 50 \ {a}. Clearly, if is nonnegative because are nonnegative. 
Therefore, H should be proportional to the Poisson kernel in f2 having mass at a, i.e.. 



H = II P(0;a,b) and F = ij,^2idP^n.a,b) for some /i G M. 

Note that is uniformly bounded from oo and 0, since H^{c^) x 1 uniformly in 6. 
To finish the proof, we need to show that fi = 1. For each 0<a^7<^t, we have 

F^{z) =^ jj ■ (1 + 0(7)) uniformly for z G [-27,27] x [0,7], 

as 5 ^ 0. Recall that D i?^(s, t) and dQ^ D L^(s) for all 6 (see ([53]) and Fig. EB)- 
Set Fq := F*^ — fiJ^^, where the function J^^ is defined in Theorem 15.41 Then Fq is 
s-holomorphic in -R^(s, t), satisfies the boundary condition (B°) on the lower boundary, 
and 

Fq{z) ^ 0(7) uniformly for z G [—27,27] x [0,7], 
smce J^^ =^ 1. Moreover, due to fCT . we have Fq{z) = 0(1) everywhere in the 
rectangle R^^{^s, ^t). Let Hq := J^{FQ{z))'^d^z, where the additive constant is chosen 
so that Hq = on the boundary L^{s). Then 

H^o = 0(a+7=^) + 05^0(1) uniformly on the boundary of Rr{2^, 7). 

Since the subharmonic function -f^olr vanishes on L^{s), Lemma [A. 31 gives 

H'oibL) < 0{6^-') ■ [0{aW) + os^o{l)] = 5 ■ 0{r'[o^ + o,^o(l)]+f ), 
where hf^^ G F denotes the inner vertex near . On the other hand, 
^o(&fnt) ^ ^\F'o{h')\' = 6\{l-^,):F\h')\' X 5I1-/.P 

which doesn't depend on a and 7. Successively passing to the limit as 5 — > 0, a — )■ 
and 7 — )■ 0, we obtain /i = l. Thus, F^ =^ \/2idP(^fi-a,b) as 5 = 5^ — 0. □ 
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(B) (D) 

Figure 9. (A) Discrete quadrilateral fi^ with four marked boundary 
points and Dobrushin-type boundary conditions. If we draw the addi- 
tional external edge [c^d^] (or [a ^d^]), there is one interface 7"^, going 
from (or rf*^, respectively) to . Since 7'' has black vertices to the 
left and white vertices to the right, its winding on all boundary arcs is 
defined uniquely due to topological reasons. (B) Besides loops, there are 
two interfaces: either and ^ d^, or ^ d^ and ^ c^. 

We denote the probabilities of these events by and Q^, respectively. 

(C) The external edge [c^d^] changes the probabilities: there is one ad- 
ditional loop (additional factor a/2), if is connected directly with a^. 

(D) The external edge [a^d^] changes the probabilities differently. 

6. 4-POINT CROSSING PROBABILITY FOR THE FK-ISING MODEL 

Let fi^ C <0 be a discrete quadrilateral, i.e. simply-connected discrete domain com- 
posed of inner rhombi z G Int and boundary half-rhombi ( G (9fi^ , with four marked 
boundary points a\ b\ c\ d^ and alternating Dobrushin-type boundary conditions (see 
Fig.lH]): i9fi^ consists of two "white" arcs a^^b^^, (Ai^i, and two "black" arcs 6^c^, (i^a^. In 
the random cluster language it means that the four arcs are wired/free/wired/free, and 
in the loop representation this creates two interfaces that end at the four marked points 
and can connect in two possible ways. As in Section 12.11 we assume that 6^ — 6^ = iS. 

Due to Dobrushin-type boundary conditions, each configuration (besides loops) 
contains two interfaces, either connecting to b^ and to d^ , or vice versa. Let 

P-^ = a\ b\ /) := ¥{a^ ^ b^; ^ d^) 

denote the probability of the first event, and Q*^ = 1 — P*^. 
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Theorem 6.1. For all r,R,t > there exists e{6) = e{6,r,R,t) such that 

if B{0,r) C C B{0,R) and either both uj{0;n^; a%^),uj{0;n^; c^d^) 
or both ijj{0;fl^;b^c^),(x!{0;fl^;d^a^) are ^t (i.e., quadrilateral has no 
neighboring small arcs), then 

\F\nl;a\b\c\d^) - p{n^;a\b\c\d^)\ ^ e{6) ^ as 6 ^ 

( uniformly with respect to the shape of and (}^), where p depends only on the con- 
formal modulus of the quadrilateral {Q^; a^, b^, c^, d^). In particular, for u G [0, 1], 

p(H;0,l-M,l,oo) - v^T^ 



-u 



Remark 6.2. This formula is a special case of a hypergeometric formula for crossings 
in a general FK model. In the Ising case it becomes algebraic and furthermore can be 
rewritten in several ways. It has an especially simple form for the crossing probabilities 



p{<P) ■.= pij])--e''^,e-'^,e'^,-e-'^) and p(f - 0) = 1 - 
in the unit disc D (clearly, the cross-ratio u is equal to svo? (p). Namely, 

p{(f)) _ sin f 



for [0,f]. 



2 ^i^HI 2' 

Curiously, this macroscopic formula formally coincides with the relative weights cor- 
responding to two different possibilities of crossings inside microscopic rhombi (see 
Fig. [HA ) in the critical FK-Ising model on isoradial graphs. 

Proof. We start with adding to our picture the "external" edge connecting and d^ 
(see Fig. [9]). Then, exactly as in Sect. 12. (12. 3 p and (12. 5p allow us to define the 
s-holomorphic in function F^^^j : fi^ — )• C such that 



on (9^2^, where 



FUO II (^(0)-^ (6.1) 



r{Q = u,{Q-u,{Q, Ce{b'c')u{d'a'), ui,2(C)er, (^-2) 

is the "discrete tangent vector" to dQ^^ oriented from /c^ to ¥ /d^ (see Fig. [9]). 
Note that F^^^^{b^) = (26)'^ and F^^^{a^) = (25)"^ ■ e-t-i'^ding(fe^-^a^)^ ^iut 

F'(d') = (26)-^ ■ g-^winding(6^^(c^^d^)) 

(and similarly for F^^^^ic^)., since the interface passes through d^ if and only if b^ is 

connected with c^, see Fig. IH]). 

Similarly, we can add an external edge \a^d^\ and construct another s-holomorphic 
in function F^^^j satisfying the same boundary conditions (16. ip . Arguing in the same 
way, we deduce that F^^^^^{b^) = {25)"^, F^^^^{c^) = (25)"^ ■ e-i«i^dmg(6^-^c^)^ ^^^^ 

F'Jd') = (26)--^ ^1 g-twinding(6^^(a*^d^)) 
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(and similarly for F^^^-^{a^)). Note that 
Let 

P^(v/2P^+Q^) ■ + Q^(P^ + v/2Q^) ■ F^^,^ 
P^(y2P'5 + Q^) + Q'^(PHy2Q'5) 
Then, also satisfies boundary conditions (16.11) . (16. 2p and, in view of (16. 3p . 

where 

^^s _ P'^(V2P^+Q^) + Q'^P'^ ^5 _ P'^Q'^ + Q^(P^+V2Q'^) 



(6.4) 



P^(v^PHQ^) + Q^(PHv^Q^) ' P^(V2PHQ^) + Q^(PH V2Q 

Since F^ is s-holomorphic and satisfies the boundary conditions (16. ip . (16. 2p . we 
can define := j^{F^{z)Yd^z and use the "boundary modification trick" (see Sec- 
tion |XS]). Then, (16. 4p implies 

i/^ = on {a^h%, H^, =0 on (a^b^), 

H^ = l on [h^c^), H^, = 1 on {b^c%,, (6.5) 

i/^ = x-^ on {c^d^) U (ci'^a'^)f , i^r* = on (c^rf'')f. U (d^a^). 



where 



x-^ = (A'^)2 = 1 - {&)^ 



{t^Y + y/2t^ 
{t^Y + ^/2t^ + l, 



p(5 p<5 



Suppose that \P\VlI) -p{VL^)\ ^ £o > for some domains Vt^ = (fi^; a\ h\ c\ d^) 
with (5 — )■ 0. Passing to a subsequence (exactly as in Section H]), we may assume that 

{Q.^; a\ h\ c\ d^) ^(fi; a, 6, c, rf), x^ ^ x G [0, 1], 



^ iJ, and F^ ^ F = V2idH, 

uniformly on compact subsets, for some harmonic function H : Q ^ M.. It follows from 
our assumptions that 5(0, r) C Q C B{0,R) and either a b, c ^ d ot b ^ c, d ^ a. 

We begin with the main case, when the limiting quadrilateral [Q; a, b, c, d) is non- 
degenerate and < X < 1. As in Section HJ we see that 

= on (ab), H = 1 on (be) and H = x on (cd) U (da). (6.7) 

Consider the conformal mapping $ from Q onto the slit strip [M x (0; 1)] \ (zx — oo; ix] 
such that a is mapped to "lower" — oo, b to +oo and c to "upper" — cxd (note that such 
a mapping is uniquely defined). Then, the imaginary part of $ is harmonic and has 
the same boundary values as H, so we conclude that if = Im$. We prove that 

d is mapped exactly to the tip %k. 
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Figure 10. (A) If d is mapped on the lower bank of the cut, then 
H < >c somewhere near [cd). The contour C = [p; q] U [q; r] U [r; s] is 
chosen so that H < x on C, dist(g; dQ) = \q—p\ and dist(r; dQ) = \r—s\. 
(B) Since = on C^, we have d^V^ ^ everywhere on C^. Moreover, 
d^V^ ^ const < on (q^r^) and d^V^ = 0([dist(u; dn^)f-'^) near p\ 
(here denotes the discrete derivative in the outer normal direction). 



Then, x can be uniquely determined from the conformal modulus of {Q; a, b, c, d). 

Suppose that above is not the case, and d is mapped, say, on the lower bank of 
the cut. It means that H < x near some (close to d) part of the boundary arc {cd). 
Then, there exists a (small) contour C = [p; q] U [q; r] U [r; s] C Q such that H < k 
everywhere on C, 

s,p G {cd) C dVL, dist(g;9i7) = \q—p\ and dist(r;9i7) = |r — s| 

(see Fig. fTOl) . Denote hy D dVl the part of VL lying inside C . For technical purposes, we 
also fix some intermediate points s, p G {sp) C {cd) and g, r G [g; r] (see Fig. [10]). For 
sufficiently small 5, we can find discrete approximations ,p^ ,p^ G {c^d^) C dQ~ 
and q^,q^,r^,r^ G fif to these points such that the contour [p^;q^] U [g''; r^] U [r"^; s^] 
approximates C. Denote by C the part of lying inside C and by C fif 
the set of all "black" vertices lying in and their neighbors. 

Let ifp := — and := Ci;^( ■ ; (s^p^); D^). Since ifp is subharmonic and 
l^'' ^ is harmonic, the discrete Green's formula gives 

[{H^{u)-H^{u,nt))V\u) - {V'{u)-V\u;nt))H^{u)]taneuu,^, > 0. 

Note that = on {s^p\ and V\u) = on := \ (sV)p. Thus, 

V'{ui^t)H^{u)tane^u,^^ ^ ^ (^^(«i,t) -i^r(«))V^'(^) tan^^„„,^, ^ 0, (6.8) 
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since H^^Ui^t) ^H^{u) everywhere on (s^p^) due to the boundary conditions fl6.1l) . fl6.2p . 
On the other hand, on most of C^, we have < (since 

on C by assumption), and V^{ui^t) ^0 everywhere on C^, which gives a contradiction. 
Unfortunately, we cannot immediately claim that < near the boundary, so one 
needs to prove that the neighborhoods of and cannot produce an error sufficient 
to compensate this difference of signs. 

More accurately, it follows from the uniform convergence =^ Co'( ■ ; (sp); D) > 
on compacts inside D and Lemma IA.3I that V^{uint) ^ const (D) ■ 5 everywhere on 
{(fr^) C C^, so, for small enough 6, 

J2 V\u,^t)Hl{u) ianeuu,^, ^ - const(D, H) < 0. (6.9) 

Thus, it is sufficient to prove that the neighborhoods of p^ and cannot compen- 
sate this negative amount which is independent of 5. Let u G {p^q^) C and 
/i = dist(u; p*^) = dist(u; dVL^) be small. Denote by i?^ the discretization of the /i x i/x 
rectangle near u (see Fig. [TOl) . Due to Lemma [A. 3| we have 

uj\u,^,-dR'^\{p'q')-Rl) = 0{5ii-'). 

Furthermore, for each v G 9/2^ \ (p^q^), Lemma [A. 21 gives 

oo^iv; (s^ify, D^) = 0(/) uniformly on 9i2^. 

Hence, 

V\u,^t) = a;'(«int; {s'j^y, D') ^ const(D) ■ 6fi~^'~^\ 
Recalling that = 0(1) by definition and summing, for any G {p^q^) C 
sufficiently close to p^ we obtain 

V\u,^,)Hl{u) tan0,,,„, ^ const(D) ■ 5 ■ (dist(M; p^))-(i-^) 

^ const(D) ■ (dist(p^;/))'^ (6.10) 

(uniformly with respect to 5). The same estimate holds near . Taking into account 
H^{u) < which holds true (if 6 is small enough) for all m G C"^ lying /x-away from p^, 
we deduce from dSTHD and fl^TTU]) that < const (D, if) ^ const (D) ■ /i^ for 

any /i > (and sufficiently small 6 ^ S{fi)), arriving at a contradiction. 

All "degenerate" cases can be dealt with in the same way: 

• if the quadrilateral {Q; a, b, c, d) is non-degenerate, then 

— if X = 1, then H < x near some part of (erf), which is impossible; 

— if (i is mapped onto the upper bank of the slit or x = 0, then H > x near 
some part of (ad), which leads to a contradiction via the same arguments 
as above; 

• if 6 = c (and so c 7^ rf), then x = since otherwise H < x everywhere near 
{cd) due to boundary conditions (16. 7p : 

• ii d = a (and so c ^ d, a then again x = since otherwise (16.71) implies 
H < K near some part of {cd) = (ca) close to a; 

• finally, a = b or c = d lead to x = 1 (otherwise H > x near some part of {da)). 
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Thus, d is mapped to the tip and so x = x{VL\ a, 6, c, d) is uniquely determined by the 
conformal modulus of the quadrilateral (x is either or 1 in degenerate cases). Recall 
that x^ = i{P^), where the bijection ^ : [0, 1] [0, 1] is given by (16. 6p . Let 

p{VL] a, 6, c, (i) := ,^~"^(x(f2; a, 6, c, d)). 

Then (since x(-) is Caratheodory stable) both P^(r2^) = ,^~^(x^) and p(f2'') tend to 
p{yt) = as 5 = 4 ^ 0, which contradicts to |P^(^^^) - pi^^)\ ^ eo > 0. 

Finally, the simple calculation for the half-plane (H; 0, 1 — u, 1, oo) gives 

H{z) ^ u + aTg[z — (1 — n)] + u arg z + (l — u) a.Tg[z — 1]), z eM. 

Hence, x(EI; 0,1 — u, 1, oo) = u and p = ^^^{u) which coincides with (II. ip . □ 

Remark 6.3. In fact, above we have shown that the "{t{z))^^ " boundary condition 
/i3.16\) reformulated in the form df^H^ ^ remains valid in the limit as 5 — )■ 0. Namely, 
let a sequence of discrete domains fl^ converge to some limiting Q in the 
Caratheodory topology, while s-holomorphic functions defined on 
satisfy l{3.16\) on arcs (s^p^) converging to some boundary arc (sp) C dQ. 
Let their integrals = Im J^{F^{z))'^d^z are defined so that they are 
uniformly bounded near (s^p^) and their (constant) values on (s^p^) 
tends to some x as 5 — )■ 0. Then, if converge to some (harmonic) 
function H inside Q, one has dnH^ ^ on (sp) in the following sense: 
there is no point ( G {sp) such that H < x in a neighborhood of (. 
The proof mimics the corresponding part of the proof of Theorem \6.1\ Since the bound- 
ary conditions i\3.16\) are typical for holomorphic observables in the critical Ising model, 
this statement eventually can be applied to all such observables. 

A. Appendix 

A.l. Estimates of the discrete harmonic measure. Here we formulate uniform 
estimates for the discrete harmonic measure on isoradial graphs which were used above. 

Lemma A.l (exit probabilities in the disc). Let uqET , r ^ 6 and a G dBf^[uQ,r) . 
Then, 

u\uo;{ay,Bl{uo,r))^6/r. 

Proof. See |Buc08j (or |CS08j Proposition A.l). The proof is based on the asymptotics 
flA.ip of the free Green's function. □ 

Lemma A. 2 (weak Beur ling-type estimate). There exists an absolute constant 
f3 > such that for any simply connected discrete domain fip, point u G Intfip and 
some part of the boundary E C dVL^ have 

"dist(M;anf)^'^ 



oj\u]E-Vt''^) ^ const' 



distfj^ (m;E) 



Here dist^a denotes the distance inside fip. 

Proof. See |CS08] Proposition 2.11. The proof is based on the uniform bound of the 
probability that the random walk on P crosses the annulus without making the full 
turn inside. □ 
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Finally, let i?p(s,t) C F denote the discretization of the open rectangle 

i?(s,t) = X (0;t) C C, s,t>0; 

G 9EIp be the boundary vertex closest to 0; and Lf.{s), U^{s, t), V^{s, t) be the lower, 
upper and vertical parts of the boundary (9i?f respectively. 

Lemma A. 3 (exit probabilities in the rectangle). Let s ^ 2t and t ^ 25. Then, 
for any = x+iy G -Rf (s,/!:), one has 



t+25 



y 



+ {y+25){t+25-y) 



and 



uj\v';V^{s,t);R^{s,t)) ^ 



t+25 s2 
x^ + {y+25){t+25-y) 



Proof. See |CS08j Lemma 3.17. The claim easily follows from the maximum principle 
for discrete harmonic functions. □ 

A. 2. Lipschitzness of discrete harmonic and discrete holomorphic functions. 

Proposition A. 4 (discrete Harnack Lemma). Let mq G F and H : B^(uo, R) — )• M 
be a nonnegative discrete harmonic function. Then, 
(i) for any Ui,U2 G Bf.{uo, r) C Int Bf,{uo, R), 



exp 

(a) for any Ui ~ Uq, 



— const ■ 



R-r 



H(U2) 

Hiui) 



const 



R-r 



\H{ui) - H{uo) I ^ const ■5H{uq)/R . 



Proof. See |Biic08j (or |CS08j Proposition 2.7). The proof is based on the asymptotics 
(lA.ip of the free Green's function. □ 

Corollary A. 5 (Lipschitzness of harmonic functions). Let H he discrete har- 
monic in B^{uo, R) andui,U2 G B^(uo,r) C Int Bf^(uo, R) . Then 

M\u2 — Ui\ 



\H{u2) — H{ui)\ ^ const 



R-r 



where M 



max \H[u) 



In order to formulate the similar result for discrete holomorphic functions we need 
some preliminary definitions. Let F be defined on some part of <). Taking the real and 
imaginary parts of d^F (see (13. 2p ). it is easy to see that F is holomorphic if and only 
if both functions 



[BF]i, 



Proj 



F(z) ; Ui(z)-U2(z) 



[WF]{i 



Proj 



F{z) ; Wi{z)-W2{z) 



are holomorphic, where Ui^2{z) G F and Wi^2{z) G F* are the black and white neighbors 
of z e <}, respectively (note that F = BF + WF). 

Let fif be a bounded simply connected discrete domain. For a function G defined 
on both "boundary contours" B,W (see Fig. [11]), we set 

„5 n—l m—1 

i G{Qd^C ■= {^{Us+l + Us)) (Us+l-Us) + {l{Ws+l+Ws)) (Ws+l-Ws). 



s=0 
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Us+l 




Let f2p be a bounded simply con- 
nected discrete domain. We denote 
hj B = uoUi..Un, Us G r, its closed 
polyline boundary enumerated in the 
counterclockwise order, and hj W = 
WoWi..Wm, Ws G r*, the closed poly- 
line path passing through the centers 
of all faces touching B enumerated in 
the counterclockwise order. In order 
to write down the discrete Cauchy for- 
mula (see Lemma lA.6|) . one needs to 
"integrate" over both B and W. 



Figure 11. Discrete Cauchy formula, notations (see Lemma [A. 6^ . 



Lemma A. 6 (Cauchy formula), (i) There exists a function (discrete Cauchy kernel) 
K{ ■ ; ■ ) : A X ^ — )• C, K{v, z) = 0{\v — z\~^), such that for any discrete holomorphic 
function F : — t- C and zq G fi^ \{B U W) the following holds true: 

F{z,) = ^.f Kiv{C);zo)F{C)d'C, 
4^ Jbuw 

where ( ~ v{C) EW,ifCeB, and C ~ ^^(C) eB,ifCeW (see Fig. MB), 
(a) Moreover, if F = BF , then 

F{zo) = FTOi—i ; ui{zo)-U2{zo) + O -— , 

where d = dist{zQ,W) , M = ma^Xz^Buw \F{z)\ and L is the length of B U W. The 
similar formula (with wi{zo) —W2{zo) instead ofui{zQ)—U2{zo)) holds true, if F = WF. 

Proof. See |CS08j Proposition 2.22 and Corollary 2.23. The proof is based on the 
discrete integration by parts and asymptotics of the discrete Cauchy kernel K{ ■ ; ■ ) 
proved by R.Kenyon in |Ken02] . □ 

Corollary A. 7 (Lipschitzness of holomorphic functions). Let F : B^(zo, R) — )• C 
be discrete holomorphic. Then there exist A, B E C such that 

F{z) = Proj[A; Ui{z) - U2{z)] + Proj[5 ; Wi{z) - W2{z)] + 0{Mr/{R-r)), 

where M = max^g^a^^^j^^^) l-^(-2)|; for any z such that \z — zo\ ^ r < R. 

Proof. Namely, 

A. 3. Estimates of the discrete Green's function. Here we prove two technical 
lemmas which were used in Section 1X51 Recall that the Green's function Gqs{-]u) : 
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fip — 7- M, n G fip C r, is the (unique) discrete harmonic in fip \ {u} function such that 
Gqs^ = on the boundary dfl^ f^ri'^) ' = 1- Clearly, 

where Gr is the free Green's function and is the unique discrete harmonic in f2f 

function that coincides with Gr on the boundary dQf^. It is known that Gr satisfies 
( |Ken02j . see also |CS08j Theorem 2.5) the asymptotics 

Griv;u) = ^ log \v-u\+o( ^ \ v ^ u. (A.l) 



Lemma A. 8. Let = i?p(2;o,r) C F, r ^ const -5, be the discrete disc, u G -Bf he 
such that Im — ^ol ^ \f 
Green's function. Then 



such that Im — ^ol ^ \f ond G = Gq6{-;u) : Bf^ M be the corresponding discrete 



\\G\\,^B^^=J2f^ii^)\G{v)\^ const -r'. 



Proof. It immediately follows from (lA.ip that 

(27r)-i log(ir) + 0{6yr') ^ G^. (■; u) <: {2n)-' log(|r) + 0(5V^') (A.2) 
on the boundary dBf^, and so inside Bf^. For G -Bp such that 1^; — -ul ^ |r, this gives 

G{v) = Gr{v] u) - G*^s {v; u) ^ -{2n)-^ log 2 + ©(^Vr^) ^ - const . 
Thus, ^ > const -r^. □ 

Lemma A. 9. Let B^ = Bf^{zo,r) C F, r ^ const -S, be the discrete disc, u G -Bp and 
G = Gqs (■; u) : -Bf — i- R 6e the corresponding discrete Green's function. Then 

ll^'^?ili,B^(.o,fO = E ^^%{m^'G]{z)\ ^ const -r. 

zeBl{zo,lr) 

Proof. Let |m — 2;o| ^ |r. It easily follows from f lA.l|) and Corollary IA.5I applied in the 
disc -B^(2;, -ul) that |[9''Gr](-2; w)| ^ const ■|z—'u|^-'^. Therefore, 

ll^'^G'r||i,ij*(^o,|r) ^ ll^''G'r||i,B^(„,2r) < const -r. 
Furthermore, double-sided bound (lA.2p and Corollary IA.5I imply 

[d^G*gs\{z) = d^[G*^s^ - (27r)-^ logr](z) = 0(const -r"^), \z-zo\ ^ |r. 

Thus, 

ll^^'^B* II 1,5^(^0, |r) ^ const -r. 
Otherwise, let j-u — zqI > f^- We have 

G*^,(-;n)^(27r)-Mog(2r) + 0(l) 

on dB^, and so on the boundary of the smaller disc B^{zo, ^r) which still contains 
B^{zq, |r). At the same time, 

Gr(-; u) ^ (27r)-i log + 0(1) on dBf.{zo, ^r). 
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Thus, 



^ G = Gr{-; u) - (■; u) ^ - const 



on the boundary, and so inside Bf^{zo, ^r). Due to Corollary IA.51 this gives 

[d^G]{z) = 0(const-r-^), for \z-zo\ ^ fr. 
Hence, \\d^G\\-^^Bi{zo,h') ^ const -r. 



□ 
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